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ABSTRACT 


— The  unimoment  method  is  applied  to  solve  the  electromagnetic 
scattering  by  a  buried  dielectric  finite  cylinder  simulating  a  land 
mine.  Computed  results  are  reported  at  frequencies  from  400  MHz  to 
1400  MHz  at  100  MHz  intervals.  The  dielectric  constants  of  the 
ground  are  considered  to  be  dispersive  which  simulates  soil  with 
5%,  10%  and  20%  water  content.  Results  are  computed  for  the  scat¬ 
tered  electric  and  magnetic  fields  which  are  presented  in  terms  of 
the  cylindrical  components  of  "E  and  ngH  at  a  distance  of  1"  to 
4"  above  the  ground  at  1"  intervals.  The  numerical  results  are 
computed  along  the  positive  x-axis  for  each  azimuthal  mode.  The 
fields  at  points  on  the  positive  x-axis  may  be  obtained  by  summing 
the  modal  fields  directly.  Fields  at  points  other  than  the  positive 
x-axis  may  be  obtained  by  summing  the  modal  fields  mutiplied  by  the 
proper  azimuthal  function  Sample  results  are  given  in  the 

report  and  the  complete  data  are  stored  on  magnetic  tape.  This 
report  includes  documentation  for  the  tapes. 


I.  INTRODUCTION 


Land  mines  made  of  plastic  or  other  dielectric  material  have  caused 
many  casualties  during  recent  conflicts.  They  are  continuously  causing 
more  casualties  many  years  after  the  conflicts,  because  these  mines 
are  difficult  to  detect  and  complete  clearing  of  such  mines  is  almost 
impossible.  One  of  the  reasons  that  an  effective  detection  system  for 
dielectric  land  mines  has  not  yet  been  developed  is  the  lack  of  theoretical 
data  for  the  electromagnetic  scattering  by  buried  dielectric  bodies. 

In  order  to  generate  data  usable  for  the  design  of  mine  detection 
systems,  this  investigation  has  successfully  applied  tne  unimoment 
method  [1]  to  compute  scattering  by  a  buried  dielectric  finite  cylinder 
which  simulates  a  dielectric  land  mine.  The  computation  uses  the 
Finite  Element  Method  (FEM)  to  treat  the  boundary  conditions  of  the 
mine.  The  solution  of  FEM  is  terminated  at  a  mathematical  sphere 
by  a  set  of  analytical  expansions  which  satisfy  the  continuity  conditions 
of  the  tangential  components  of  E  and  H  fields  on  the  air-ground 
interface.  Owing  to  the  versatility  of  the  FEM  in  matching  the  mine 
surface  and  the  fast  convergence  of  the  analytical  expansions,  we  are 
able  to  compute  the  scattering  by  the  buried  dielectric  land  mines 
efficiently. 

This  report  provides  the  computational  results  for  eleven  different 
frequencies:  400  HMz  to  1400  MHz  at  100  MHz  interval.  Five  different 
incident  angles,  each  with  two  different  polarizations,  are  considered 
for  the  incident  plane  wave.  Field  components  are  calculated  on  the 
horizontal  planes  at  four  different  altitude  levels,  i.e.,  l"-4"  above 
ground  at  1"  interval.  The  dielectric  constants  are  assumed  to  be 
dispersive  using  the  formulas  offered  by  Von  Hippel,  simulating  silt 
loam  with  5%,  10%  and  20%  water. 

Results  are  computed  for  the  p,  <J>,  z  components  of  both 
scattered  electric  and  magnetic  fields  for  each  azimuthal  mode  on 
the  positive  x-axis  at  1  cm  intervals  up  to  15  ans.  The  total  scattered 
fields  on  the  positive  x-axis  can  be  obtained  by  directly  summing  the 
modal  components.  The  scattered  fields  at  locations  other  than  the 
positive  x-axis  may  be  obtained  by  summing  the  modal  fields  multiplied 


by  the  corresponding  ejm<^  factor.  Therefore,  we  are  able  to  obtain 
2-dimensional  data  from  a  one  dimensional  data  set,  which  results  in 
savings  in  data  handling. 

Due  to  the  massive  amount  of  data  involved,  it  is  not  feasible  to 
present  all  of  them  graphically.  The  data  is  stored  in  magnetic  tapes 
This  report  also  contains  documentation  for  the  tapes  so  that  a 
specific  data  set  can  be  conveniently  fetched. 

Examples  are  given  of  how  to  use  the  data  to  compute  the  total 
and  scattered  electric  and  magnetic  fields.  Typical  data  sets  of 
azimuthal  coefficients  are  given  in  Chapter  VIII.  Appendix  I  gives 
examples  of  using  the  data  to  calculate  the  scattered  and  total  fields 
Appendix  II  gives  a  listing  of  program  PCOEF,  which  is  to  be  used  to 
read  selected  portions  of  the  data  from  the  magnetic  tapes.  Appendix 
III  gives  a  table  of  the  various  symbols  used  in  this  report  and  their 
meani ngs . 
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II.  THE  unimoment  method 


The  unimoment  method  was  first  published  in  1974  [1]  to 
generalize  the  moment  method  for  the  radiation  and  scattering  problems 
involving  material  bodies.  It  was  then  applied  to  the  inhomogeneously 
loaded  biconical  antennas  [2],  the  scattering  by  two  dimensional 
dielec'-n'c  cylinders  [3],  and  by  axially  symmetric  dielectrical  bodies 
[4,  5].  Recent  discovery  of  the  generalizations  of  Sommerfeld  integrals 
and  a  new  type  of  field  expansions  in  two  medium  half  spaces  [6]  have  made 
possible  the  extension  of  the  method  to  the  electromagnetic  scattering 
by  buried  dielectric  land  mines  [7]. 

The  basic  idea  of  the  unimoment  method  is  to  combine  the  Finite 
Element  Method  (FEM)  with  the  analytic  solutions  so  that  the  maximum 
advantages  are  obtained  from  both.  The  analytic  solutions  are  usually 
expressed  in  terms  of  truncated  series  expansions.  Although  analytic 
expansions  usually  represent  the  fields  of  a  vast  region  in  space, 
they  are  only  applicable  to  simply  shaped  objects  such  as  spheres.  The 
FEM  has  been  noted  for  its  flexibility  in  fitting  boudaries  with 
general  shapes.  For  scattering  problems  which  involve  an  infinite  - 
space,  it  is  impractical  to  use  the  FEM  for  the  entire  space.  To 
minimize  the  number  of  unknowns,  the  FEM  is  terminated  at  two  spheres; 
one  encloses  the  entire  mine,  and  the  other  is  imbedded  in  the  mine. 
Spherical  vector  wave  expansions  represent  the  fields  in  the  interior 
of  the  sphere  inside  the  mine.  In  the  exterior  of  the  outer  sphere, 
complications  of  the  analytic  expansions  occur  due  to  the  presence  of 
the  air-ground  interface.  The  computation  in  this  report  has  employed 
a  special  type  of  vector  wave  expansion  incorporated  with  the  generalized 
Sommerfeld  integrals  to  satisfy  the  air-ground  boundary  conditions. 

The  coefficients  of  the  analytical  expansions  are  obtained  by  enforcing 
the  continuity  conditons  on  the  spheres.  The  scattered  fields  can  then 
be  calculated  from  the  coefficients  and  the  modal  fields. 

We  shall  now  briefly  describe  the  procedures  of  the  unimoment 
method  in  the  following  steps: 


(1)  Draw  a  mathematical  sphere  S-|  to  enclose  the  entire  land  mine 
as  shown  in  Figure  1. 

The  sphere  mathematically  separates  the  space  into  two  parts. 

The  interior  of  the  sphere  {Region  III)  contains  the  land  mine.  The 
exterior  of  the  sphere  (Regions  I  and  II)  involves  the  lossy  ground 
and  the  air.  For  efficiency  and  convenience,  the  sphere  should  be  as 
small  as  possible  to  reduce  the  FEM  calculation,  and  large  enough  to 
enclose  the  entire  mine. 

(2)  Draw  a  supplemental  mathematical  sphere  S2  in  the  interior  of 
the  land  mine  as  shown  in  Figure  2. 

The  supplemental  sphere  is  drawn  to  reduce  the  region  of  FEM 
further.  The  sphere  $2  mathematically  separates  the  region  inside 
S.j  into  two  parts.  Region  IV  is  the  interior  of  S2  and  Region  III 
is  the  area  between  spheres  and  S2.  For  efficiency  and  convenience 
S2  should  be  as  large  as  possible  but  small  enough  to  be  completely 
inside  the  mine. 

(3)  Solve  the  Maxwell's  equations  for  Region  III  between  S-j  and  S2 
using  FEM. 

The  original  Maxwell's  equations  has  six  vector  components  in  3- 
dimensions.  Because  of  axial  symmetry  of  the  geometry,  it  is  possible 
to  reduce  the  six  unknowns  to  two  coupled  azimuthal  potentials 
=  kopE^  and  ^  =  k0pH^  for  each  azimuthal  Fourier  series  mode 

with  e^  variations.  All  other  components  of  ¥  and  H  fields 
can  be  directly  obtained  from  these  potentials.  The  potentials 
satisfy  two  coupled  differential  equations  in  Region  III  including 
the  boundary  of  the  mine. 

By  using  the  Fourier  series  of  e,',m<^)  and  the  coupled  azimuthal 
potentials  we  are  able  to  reduce  the  original  three  dimensional  problem 
into  many  much  smaller  two  dimensional  problems.  The  amount  of  computer 
time  can  thus  be  greatly  reduced. 

The  differential  equations  are  then  changed  to  their  corresponding 
variational  integrals.  The  functions  of  the  potentials  which 
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Figure  2.  Minimization  of  FEM  region  III  by  a 
supplemental  inner  sphere  S9 


L 


render  stationary  the  variational  integral  are  also  the  solutions  of  the 
differential  equations.  The  basic  idea  of  FEM  is  to  approximate  the 
potentials  by  piecewise  polynomials  so  that  the  variational  integrals 
can  be  evaluated  analytically. 

In  the  FEM,  we  first  divide  the  entire  Region  III  into  many 
small  trianglar  elements  such  as  those  shown  in  Figure  3.  The 
triangles  conform  with  the  surfaces  of  the  mine,  the  inner  sphere, 
and  the  outer  sphere.  The  elements  inside  the  mine  have  relative 
dielectric  constant  em  =  2.89  and  those  outside  the  mine  have  ee. 

The  subscripts  m  and  e  denote  mine  and  earth  respectively. 

The  potentials  in  each  triangle  are  then  represented  by  inter¬ 
polating  polynomials  passing  through  each  nodal  value.  The  variational 
integral  is  obtained  by  summing  up  the  element  integrations  over 
all  elements  in  Region  III.  The  variational  formulation  is  stationary 
if  its  differentiation  with  respect  to  nodal  values  is  zero.  This 
leads  to  a  set  of  linear  equations  which  can  be  used  to  solve  for 
the  fields  in  Region  III.  The  matrix  involved  in  these  equations  is 
a  banded  sparse  matrix  which  can  be  solved  efficiently  by  a  special 
sparse  matrix  algorithm. 

Since  the  boundary  values  on  the  spheres  are  not  specified,  the 
solutions  of  Region  III  are  in  general  not  unique.  In  fact  they  result 
in  a  set  of  linearly  independent  solutions. 

(4)  Expand  the  field  in  Region  IV  inside  the  sphere  Sg  by  using  the 
conventional  spherical  vector  waves. 

Ihe  field  in  Region  IV  is  expanded  analytically  by  the  conventional 
spherical  vector  waves.  Since  the  origin  of  the  coordinate  is  included 
in  this  region,  only  the  spherical  Bessel  Functions  of  the  first  kind 
jn(kr)  is  required  as  the  radial  function? 

(5)  Represent  the  field  in  Region  I  and  II  outside  the  sphere  S-j  by 
mutipole  expansions  which  satisfy  the  air-ground  boundary  conditions. 
The  fields  in  Regions  I  and  II  are  decomposed  into  incident  and 

scattered  fields.  The  scattered  fields  in  Region  II  (the  L'^th)  are 
further  decomposed  into  direct  and  secondary  waves. 
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The  direct  field  is  represented  in  terms  of  multipole  expansions 
as  if  the  scatterer  is  in  an  infinite  medium.  The  secondary  fields  are 
represented  in  generalized  Sommerfeld's  integrals.  The  sum  of  the 
direct  and  secondary  fields  satisfies  the  air-ground  interface  boundary 
conditions.  The  rectilinear  spherical  vector  waves,  in  which  the  rec¬ 
tangular  vectors  and  the  spherical  harmonics  are  combined,  form  the  vector 
potentials.  This  type  of  combination  was  introduced  to  reduce  the  com¬ 
plexities  in  enforcing  both  the  spherical  and  planar  boundary  conditions 
on  the  mathematical  sphere  and  on  the  earth  surface.  Supplemental  terms 
using  horizontal  rotating  multipoles  have  been  discovered  and  added  to 
the  direct  field  expansions  to  speed  up  the  convergence  [6], 

The  direct  fields  are  then  transformed  into  cylindrical  harmonics 
by  using  the  Fourier-Bessel  integrals  which  replace  the  spherical-Hankel- 
Lengendre  functions.  By  enforcing  the  continuity  conditions  of  the 
tangential  E  and  H  fields  on  the  air-ground  interface,  the  secondary 
fields  in  Region  II  and  the  scattered  field  in  Region  I  are  found  to  be 
in  the  form  of  "Generalized  Sommerfeld  integrals"  [6],  The  expansions 
with  the  generalized  Sommerfeld  integrals  satisfy  the  air-grcund  boundary 
condition  term-by-term. 

The  analysis  has  thus  yielded  a  converging  multi  pole  expansion  tech¬ 
nique  for  the  fields  in  two-medium  half  spaces.  The  same  set  of  expansion 
coefficients  are  used  in  both  Regions  I  and  II. 

(6)  Solve  the  expansion  coefficients  by  enfoi..;.\g  the  continuity 
conditions  on  both  the  outer  sphere  S-j  and  the  inner  sphere 
S2. 

The  three  sets  of  solutions  in  Regions  II,  III,  and  IV  are 
solved  individually  by  the  previously  described  methods.  We  then 
couple  all  the  solutions  on  the  mathematical  spheres  S-j  and  S^ 
by  using  the  continuity  conditions  of  tangential  electric  and  magnetic 
fields.  The  incident  waves  are  the  driving  fields  for  the  calculations. 

(7)  Calculate  scattered  fields  in  Region  I  (the  air)  by  using  the 
coefficients  and  modal  fields. 

The  final  step  of  the  unimoment  method  is  to  generate  the 
required  scattered  electromagnetic  fields  by  a  buried  land  mine. 

The  fields  on  planes  parallel  to  the  ground  at  different  heights 
above  the  ground  are  calculated. 


c 
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III.  GROUND  PARAMETERS 

In  an  earlier  computation,  [8],  the  relative  dielectric  constant 
of  the  ground  was  assumed  to  be  constant  at  eg  =  9-j7  for  all  frequencies. 
This  unrealistic  assumption  has  been  corrected  in  this  report.  The  complex 
dielectric  permittivities  are  now  considered  to  be  functions  of  frequencies 
as  given  by  Falls  and  Mittleman  [10]. 

In  this  report,  we  consider  three  types  of  soil  conditions,  i.e., 
silt  loam  with  5%,  10%  and  20%  water.  The  real  part  of  er  and  the 
attenuation  parameter,  a  ,  for  the  soil  are  given  in  Table  I. 

The  conversion  from  er  and  a  to  er  and  (real  and 
imaginary  part  of  Eg)  is,  (eq.  3.38  in  ref.  [9]), 


The  real  and  imaginary  part  of  e0  which  are  used  in  the  computations 
for  this  report  are  tabulated  in  Table  II.  In  equation  (1),  tan  6 
is  called  the  loss  tangent  and  is  the  free  space  wavelength. 


TABLE  1 

DIELECTRIC  CONSTANTS  [SILT  LOAM] 


Frequency 

1  'i - — — 

5%  Water 

10%  Water 

20%  Water 

1 

[Mega  Hertz] 

a 

[db/meter] 

er 

a 

[db/meter] 

1  — 1 . . 

er 

a 

[db/meter] 

e 

r 

400 

15.5 

5.2 

27.0 

7.8 

30.0 

17.6 

500 

19.0 

5.1 

33.0 

B 

34.0 

17.0 

600 

20.2 

5.0 

35.0 

B 

37.5 

16.8 

700 

22.0 

a 

38.0 

7.7 

40.5 

16.6 

80C 

23.8 

B 

40.5 

7.8 

43.5 

-16.4 

900 

25.5 

4.8  1 

43.0 

7.8 

46.5 

16.2 

1000 

27.0 

4.7 

45.5 

7.8 

49.2 

, 

16.0 

1100 

29.2 

B 

48.0 

7.8 

52.0 

'  15.9 

1200 

31.8 

50.7  ; 

7.7 

54.9 

15.8 

1300 

34.0 

■ 

53.0 

• 

7.6 

57.5 

15.7 

1400 

36.5 

B 

56.0 

7.6 

60.2 

15.6 

n- 
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TABLE  II 

DIELECTRIC  CONSTANTS  [SILT  LOAM] 


PR  ECU  ENCY 
(  PC GA  HERTZ) 

400  .00 
500. CC 
6C0.CC 
700  .00 

ecc. cc 

900. CC . 

I ccc.co 
1  ICC. cc 
12CC.C0 
I3CC.00 
1  4  C  C  .  C  0 


FPL'  CCt  NC  Y 
(MFC*  HERTZ  ) 

4C0. CO 
5C0.C0 
600  .00 
7CC. CO 
8C0.0C 
900  .00 
1CCC.  CO 

..  .  .lic-c. co 

1  2  CC .CO 
13CC.CC 
1400  .00 


L  OS  S 

(DUm'UM 

15  .50 
19.  CO 
2  0.2C 
22  .00 
23.  ec 
2  5  .  i  C 
27  .CO 
29.20 
31  .CC 
34  .00 
36.  5 C 


LCSS 

(  CU/ME  TER) 

30.CC 
.  3A.CC 

37.50 
AO.  5 C 
43. 5C 

46 .50 
49. 2C 
52.CC 
54  .90 
57.  5 C 
6"  ■  ?0 


REAL  PAM 
or  THE 
r.FLAT  I  vi; 

PERM1  m  Vi  7  V 
5.20 
5.10 
5.  CC 
4.90 

4 . 80 
4.  EC 

4.70 

4.60 
4.  6  0 
A  .  50 
4.  60 


REAL  PART 
OF  TEE 

relative 

PERM1TT1  VI  TV 

7.80 

7.60 
7.  6  0 

7.70 
7.  EC 
7.  80 
7;  00 
7.  EC 

7.70 

7.60 
7.  6 C 


real  part 
■  Cf  The¬ 
ft  C  LA  71  VE 
PERM  ITT  I  V  IT  Y 
1  7  .60 
1  7,  CC 

16 .80 
1  6  .60 
1  6.  4C 
16  .,20 
16  .bo 
15.  SC 
15.00 
1  5.  ?C 
15.  6 C 
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S 1LT  LOAM 
TYPE  10  PER 


T  R FLUENCY 
(MEGA  HERTZ) 

ACC. CO 
5CC. CO 
600. CO 
7CC.CC 

tcc.cc 

SOO  .00 
1CC0. CO 
1 100.00 
1200  .00 
13CC.CC 
1AC0.CC 


SILT  LOAM 
TYPE  20  PER 


CENT  HAH  R 


LCSS 

(  OB  /  ME  TF  R  ) 

27  .CC 
33.  CC 
35.00 
30  .CO 
AO.  5C 
A3  .CO 
A 5 .  50 
A8.CC 
50 . 7C 
53.CC 
56  .ec 


CERT  WATER 


IM  AG  I  MARY 
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IV.  NUMERICAL  RESULTS 


Using  the  unimoment  method  described  in  Chapter  II,  this  investigation 
has  successfully  computed  the  scattering  of  electromagnetic  waves  from  a 
buried  object  which  is  electrically  very  similar  to  an  antipersonnel 
dielectric  land  mine.  The  geometrical  configuration  of  the  proposed 
computation  is  shown  in  Figure  4.  The  buried  target  is  a  finite  cylinder 
which  is  5.5  cm  high  and  11.2  cm  in  diameter.  The  dimensions  approximate 
those  of  an  antipersonnel  land  mine.  The  relative  dielectric  constant 
of  the  target  is  cr  =  2.89,  which  is  that  of  TNT.  The  relative  dielectric 
constant  of  the  ground  eg  as  functions  of  frequency  are  given  in  Tables 
I  and  II.  The  top  edge  of  the  mine  is  buried  4.1  cm  under  the  surface 
of  the  ground.  This  configuration  is  believed  to  be  close  to  the  actual 
situation. 

The  scattering  configuration  involves  plane  waves  incident  at 
0^  =  0°,  30°,  45°,  60°,  75°  with  respect  to  the  z-axis.  The  frequencies 
considered  are  at  400  MHz  to  1400  MHz  at  100  MHz  interval.  The  scattered 
electric  fields  are  computed  on  the  planes  parallel  to  the  ground  at 
z  =  1"  (2.5  cm),  2"  (5.08  cm),  3"  (7.62  cm),  and  4"  (10.16  cm). 

The  x  and  z  axes  are  defined  so  that  the  x-z  plane  is  the 
plane  of  incidence.  The  plane  of  incidence  is  defined  to  be  the  plane 
which  contains  the  propagation  vector  ¥  .  Two  different  kinds  of 
polarizations  are  considered  for  each  incident  angle.  They  are  t1 
in  the  x-z  plane  (or  H-Y  incidence)  and  ft1  in  the  x-z  plane 
(or  E-Y  incidence).  Owing  to  their  basic  differences  in  refraction 
and  scattering,  the  two  polarizations  should  result  in  quite  different 
scattered  field  patterns  except  for  the  symmetric  case  when  0^  =0°. 

In  all  the  computations  shown  in  this  report,  the  original  incident 
plane  waves  (before  reflection  and  refraction  from  the  ground)  are 
considered  to  have  an  amplitude  of  1  volt/meter,  i.e., 

it’j =  %  ifi1 1  =i 


The  computation  involves  5  incident  angles  each  with  2  polar¬ 
izations  at  11  frequencies.  For  each  case,  the  scattered  fields  are 
given  in  p,  <j>,  z  components  at  4  different  parallel  planes  above 
ground.  The  total  number  of  data  cases  to  be  presented  are  5  x  2 
x  11  x  3  x  4  =  1,320. 
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Because  of  the  vast  amount  of  data  involved,  it  is  impractical 
to  present  all  the  data  in  graphic  form.  They  are  recorded  on  magnetic 
tapes  so  that  MERADCOM  engineers  may  recall  the  data  at  their  convenience. 
The  documentation  of  the  tapes  are  in  Chapters  VI  and  VII  of  this  report. 

In  Figures  5  -  10,  the  typical  computed  results  are  shown  for 
600  MHz,  1000  MHz  and  1400  MHz. 

It  is  noticed  in  Figures  8-10,  where  the  soil  has  a  20% 
water  content,  that  the  scattered  fields  for  1000  MHz  is  higher  than 
both  600  MHz  and  1400  MHz.  This  phenomenon  reverses  the  trend  of 
previous  computations  [8],  which  indicates  the  scattered  fields  to 
be  higher  at  lower  frequencies.  Since  the  previous  computations 
were  based  on  an  unrealistic  assumption  that  the  ground  parameter 
was  independent  of  frequencies,  the  present  results  are  definitely 
more  reliable. 
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Figure  5  . 

Curves  show  that  1  inch  above  the  ground  the  magnitude  of  Es  qenerallv 
increases  with- increased  water  content.  '  P 
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Figure  6 

Curves  show  that  1  inch  above  the  ground  the  magnitude  of 
increases  with  increased  water  content.  1 
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MAGNITUDE  OF  ETAO*H-PHI  IN  VOLT/METER 


ANGLE  OF  INCIDENCE  IS  0.  WITH  H  POLARIZED  IN  THE  -Y  DIRECTION. 


1  .  4-OE— 01 


1 . 20E— 01 


1 . 00E— 01 


8.00E-02 


6.00E-02 


4.00E-02- 


2.00E-02 


0.00 


-15.0 


-10.0 


RHO  IN  CENTI-METERS 


PHI  CUT  OF  180  -  0  DEGREES  IS  USED. 

AT  1.0  INCH  ABOVE  THE  GROUND 

♦  SILT  LOAM  WITH  5.0  PER  CENT  WATER  FREQUENCY  600.0  MEGA  HERTZ. 

X  SILT  LOAM  WITH  10.0  PER  CENT  WATER  FREQUENCY  600.0  MEGA  HERTZ. 

+  SILT  LOAM  WITH  20.0  PER  CENT  WATER  FREQUENCY  600.0  MEGA  HERTZ. 

Figure  7 

Curves  show  that  1  inch  above  the  ground  the  magnitude  of  generally 
increases  with  increased  water  content  at  600  MHz,  ^ 
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Figure  3 

Curves  show  that  1  inch  above  the  ground  the  magnitude  of 


is  generally 


large;  at  1000  MHz  than  at  either  600  or  1400  •  IMz  for  silt  loam  having  a 
20%  water  content.  -19- 


MAGNITUDE  Of  E-Z  IN  VOLT/METER 


MAGNITUDE  OF  ETAO*H-PHI  IN  VOLT/METER 


V.  DATA  PRESENTATION 


The  data  for  the  scattered  fields  are  presented  in  such  a  way 

as  to  give  the  maximum  information  with  the  minimum  data  set  and 

easiest  access.  The  fields  are  decomposed  in  their  cylindrical 

components  Es,  E^,  E^  and  n  Hs,  n  H^,  n  H^.  For  a  particular 
p  9  z  o  p  o  9  0  z 

incident  angle,  polarization,  and  p,  a  field  location,  the 
fields  are  given  in  terms  of  their  "Azimuthal  Coefficients" 

For  example: 

M 

E^  =  V''  a  eJ"m^  C  =  P,  4>  or  z 

C  Z— *  m 

-M 


where  a  's  are  the  "Azimuthal  Coefficients".  Therefore,  along 
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-S  . 


the  positive  x-axis,  E  which  equals  E  is  as  shown  below 
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And,  along  the  negative  x-axis,  E  which  equals  -Ev  is  as  shown  below, 
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The  field  at  any  other  azimuthal  angles  can  be  found  accordingly. 

In  this  report,  the  "Azimuthal  Coefficients"  are  computed  for 
p  from  0  to  15  cm  at  1  cm  interval,  and  for  z  from  1-4"  at  1"  interval 
Therefore,  by  proper  summation  of  the  "Azimuthal  Coefficients"  fields 
within  circles  of  15  cm  radius  at  4  different  height  can  be  obtained. 

The  "Azimuthal  Coefficients"  are  presented  in  tapes,  the 
documentation  of  which  is  in’ the  next  two  chapters.  Samples  of  the 
information  on  the  tapes  are  given  in  Chapter  VIII. 


VI. 


Description  of  the  Program  PCOEF 
(sent  via  nine  track  tape) 
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1.  A  complete  listing  of  the  Fortran  Source  of  the  program  PCOEF 
is  in  Appendix  II. 

r 

2.  Here  at  Berkeley  the  MNF4  compiler  was  used.  The  input  consists  of  the 
file  of  coefficients  to  be  read,  TAPE  18,  and  the  instructions  by  the 
user  which  are  in  the  file  input  which  has  been  made  equivalent  to  TAPE 
10,  by  the  Program  statement  in  the  main  part  of  the  Program. 

3.  TAPE  18  should  be  positioned  at  its  beginning. 

4.  Below  is  a  typical  input  for  TAPE  10: 

500.0  05.0  30.0  00.0  2.0  2. 

ETA0*HPHI  EZEE 

1 2345678901 2345678901 2345678901 2345678901 2345678901 234567890 
1 1 11 1 1 11 11 22222222223333333333444444444455555555556 

5.  The  input,  TAPE  10,  is  read  by  subroutine  READIN.  The  reader  should 
refer  to  that  subroutine  now.  The  first  logical  record  contains  the 
desired  frequency  (in  Mega  Hertz),  the  pev'  cent  water  content  of  tne 
soil,  the  angle  of  incidence  in  degrees,  the  angle  of  polarization 

in  degrees,  the  height  in  inches,  and  the  number  of  components  desired. 
The  above  example  in  Item  4  indicates  500  MHz,  5%  water,  30°  incident 
angle,  0°  polarization,  2  inches  above  the  ground  and  2  field  components. 
Subroutine  READIN  will  then  call  subroutine  RCNAME;  it  reads  the 
names  of  the  components  in  the  order  they  are  to  be  printed.  Note 
the  names  must  be  exactly  as  listed  in  array  AC0MP,  which  is  set 
by  the  third  data  statement  in  the  main  part  of  the  program.  Again, 
for  the  sake  of  clarification,  the  names  must  be  the  same  as  the 
elements  in  array  ACOMP,  such  as  ETA0*HPHI  for  OqH^,  EZEE  for  Ez  etc., 

but  any  order  is  allowed.  Note,  the  first  character  of  each  element 
of  ACOMP  is  a  blank. 

6.  The  subroutine  READIN  uses  subroutine  EOFILE  to  be  sure  the  input 
which  is  TAPE  10  will  be  positioned  at  the  start  of  the  next  file. 

A  user  at  a  computer  center  other  than  LBL's  may  have  to  change 
EOFILE,  since  it  calls  EOF. 

7.  This  program  has  been  successfully  run  here  at  Berkeley. 
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VII. 


Description  of  the  files  Containing  Azimuthal 
Coefficients  sent  ia  the  nine  track  tapes 
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1.  The  tapes  were  written  by  using  LBL's  ENCODE.  A  record  length  of  80 
characters  and  a  blocking  factor  of  45  was  used.  The  tapes  are 
written  in  ASCII. 

2.  Each  file  containing  azimuthal  coefficients  is  for  a  mine  buried 
in  silt  loam  having  either  5%,  10%  or  20%  water  content.  The 
frequency  of  the  incident  fields  is  one  of  the  following:  400, 

500,  600,  700,  800,  900,  1000,  1100,  1200,  1300,  or  1400  Mega 
Hertz . 

3.  The  air-ground  interface  is  taken  to  be  the  xy  plane,  with  the 
mine  being  buried  and  centered  on  the  negative  z  axis.  The  incident 
planar  field's  Poynting  vector  is  taken  to  always  lie  in  the  xz 
plane.  Incident  angles  are  with  respect  to  the  vertical  axis 

(the  z  axis).  Plane  waves  whose  incident  angles  are  0,  30,  45, 

60,  and  75  degrees  were  used.  Note  that  an  incident  angle  of  0 
degrees  corresponds  to  a  planar  wave  whose  Poynting  vector  is 
in  the  negative  z  direction.  For  each  incident  angle,  two  polar¬ 
izations  are  considered.  The  polarizations  are  defined  by  a  polar¬ 
ization  angle.  The  angles  of  0  and  90  degrees  were  used.  A 
polarization  angle  of  0  degrees  means  that  the  incident  H  field 
is  polarized  in  the  -y  direction.  A  polarization  angle  of  90 
degrees  means  that  the  incident  E  field  is  polarized  in  the  +y 
direction. 

4.  To  obtain  the  total  fields,  add  the  incident  field  and  the  field 
scattered  by  the  air-ground  interface  to  that  defined  by  the  azimuthal 
coefficients.  The  magnitudes  are  normalized  so  that  the  electric 
field  of  the  incident  plane  wave  has  a  value  of  1  volt  per  meter. 

The  azimuthal  coefficients  for  the  phi,  rho,  and  zee  components  of 
the  electric  field  are  given  in  units  of  volts  per  meter.  Azimuthal 
coefficients  for  the  phi,  rho,  and  zee  components  of  the  magnetic  field 
are  given  in  n0H,  (n  =  377ft).  These  coefficients  are  also  in  units  of 
volts  per  meter.  ° 

5.  The  files  containing  the  azimuthal  coefficients,  when  printed,  can 
be  read  by  a  human  being,  although  this  would  certainly  be  a  tedious 
task.  The  results  are  grouped  by  incident  plane  wave.  For  each 
incident  plane  wave,  the  results  are  grouped  by  a  particular  height 
above  the  air-ground  interface.  For  each  particular  height  above 

the  air-ground  interface,  the  results  are  grouped  by  field  components. 
For  each  field  component,  the  results  are  grouped  by  a  particular 
value  of  rho.  For  each  value  of  rho,  the  most  negative  azimuthal 
order  appears  first. 

6.  What  now  follows  is  a  line  by  line  description  of  such  a  file. 

The  reader  may  now  wish  to  refer  to  a  listing  of  TAPE  18. 

The  first  part  of  Chapter  VIII  of  this  report  contains  a  listing  of 
the  first  part  of  one  of  the  files. 
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7.  First  appears  a  heading.  The  first  line  of  this  eight  line  heading  is 
as  shown  directly  below: 


FREQMH ,PCENT , Z INCH , AOF I , AOFP ,MM3 ,NPTU 


FREQMH  indicates  frequency  in  Mega  Hertz 

PCENT  indicates  per  cent  water  content  of  the  soil 

ZINCH  indicates  height  above  the  air-ground  interface  in  inches 

AOFI  indicates  angle  of  incidence  in  degrees 

AOFP  indicates  angle  of  polarization  in  degrees.  (As  stated  above, 

there  are  only  two  values  used.  0  corresponds  to  the  incident 

H  field  polarized  in  the  -y  direction;  90  indicates  the 
incident  E  field  is  polarized  in  the  +y  direction. 

MM3  indicates  the  total  number  of  azimuthal  coefficients.  This 

number  is  always  a  positive  odd  integer,  since  the  absolute 
value  of  the  most  negative  azimuthal  order  used  equals  the 
value  of  the  largest,  with  all  azimuthal  orders  lying  between 
these  limits  being  used  as  well. 

NPTU  indicates  the  number  of  rho  values  used  for  the  height. 

The  number  used  is  always  16. 


8.  The  second  line  of  the  heading  lists  the  values  of  the  first  line,  in 
the  order  of  the  first  line.  Referring  to  the  first  heading  of  the 
file  shown  in  Chapter  VIII,  we  have  freq.  =  1400  MHz,  20%  water,  1  inch 
above  ground,  0°  incident,  0°  polarization,  MM3  -  11,  NPTU  =16.  The 
third  line  prints  the  format  to  be  used  to  read  the  first  three  lines 
of  the  heading. 

9.  The  next  line  is  as  shown  below: 

(RHO ( I ) ,J=1 ,NPTU,1)  —  IN  F0RMAT(4(1X,E10.4)) 

RHO  indicates  the  rho  values  used  in  centimeters.  The  format  is  the 
format  used  to  read  the  next  four  lines,  which  contain  the  various 
rho  values  used.  The  rho  values  used  are  always  the  same  0  through 
15  centimeters  in  increments  of  1  centimeter. 

10.  These  first  eight  lines  then  are  what  are  referred  to  as  the  heading. 

A  heading  then  always  appears  before  the  results  for  each  particular 
height. 

11.  Following  the  heading  are  two  lines.  These  two  lines  indicate  the 
field  component  and  the  format  to  use  to  read  these  two  lines. 

Following  this  is  a  line  shown  below: 

( ( WAVEH ( IR0W , IC0L ) ,  IR0W=1 ,MM3,1) ,JC0L=1NPTU,1)  —  IN  F0RMAT( (6(1X,E10.4) ) 


11.  (cont.) 


This  indicates  that  the  two  dimensional  array  for  storing  the 

azimuthal  coefficients  are  written  out  in  groups  of  MM3(=11  in 

the  example  shown  in  Chapter  VIII)  complex  numbers,  (22  real  numbers). 

The  first  block  of  coefficients  refers  to  JC0L=1,  which  represents 
EPHI  at  p  =  0.  The  first  two  numbers  in  JC0L=1,  represents  the 
real  and  Imaginary  part  of  the  most  negative  azimuthal  order, 
i.e.,  M=-5  for  MM3=11.  The  next  two  numbers  correspond  to  M=-4, 
etc,  and  the  last  two  numbers  corresponds  to  M=5.  The  next  block 
of  numbers,  JC0L=2,  represents  EPHI  at  p  =  1  cm,  and  each  sub¬ 
sequent  blocks  represents  the  field  at  1  cm  increment  in  p. 

12.  Following  these  coefficients  are  two  lines  which  indicate  the  next 

field  component,  which  are  in  turn  followed  by  its  azimuthal  coefficients. 
This  pattern  is  then  repeated  throughout  the  tape. 

13.  Appendix  I  illustrates  the  application  of  the  data  base  to  compute 
scattered  and  total  fields  on  planes  above  the  ground. 


Typical  Portions  of  the  data  Base 
of  Azimuthal  Coefficients 


The  incident  field  in  this  report  refers  to  the  incident  plane 
wave.  The  primary  field  refers  to  the  incident  field  plus  the  reflected 
field  from  the  interface  without  the  buried  object.  The  scattered  field 
is  meant  to  be  the  total  field  less  the  primary  field.  The  azimuthal 
coefficients  are  for  the  scattered  fields.  To  obtain  the  total  fields, 
add  the  primary  field  to  that  defined  by  the  azimuthal  coefficients. 

The  magnitudes  are  normalized  so  that  the  electric  field  of  the  incident 
plane  wave  has  a  value  of  1  volt  per  meter.  The  azimuthal  coefficients 
for  the  phi,  rho,  and  zee  components  of  the  electric  field  are  given 
in  units  of  volts  per  meter.  Azimuthal  coefficients  for  the  phi,  rho, 
and  zee  components  of  the  magnetic  field  are  given  in  n  H,  (n  =  377ft). 
These  coefficients  are  also  in  units  of  volts  per  meter?  0 

The  coefficients  shown  are  from  the  beginning  of  a  file  where  the 
Incident  field's  frequency  is  1400  Mega  Hertz  and  the  soil  is  Silt  Loam 
having  a  20%  water  content.  The  near  field  is  for  a  distance  of  1  inch 
above  the  air-ground  interface.  Coefficients  in  Set  I  are  for  0°  incidence 
and  those  of  Set  II  are  for  30°  incidence.  Both  sets  of  data  are  for 
polarization  angle  of  zero  degree. 
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0.  0.  0. 

-.2829E-02  -.5920E-02  0. 

0.  0.  0. 


-.2783E-0I  -.2I68E 
0. 

0. 

.9300E-01  .2115E 

0. 

0. 

.  0. 

0. 

.  0. 

0. 

.  0. 


-.5016E-01  -.1198E 


•921IE-01  -.9058E 


1789E-01  0. 


0.  3. 

3.  0. 

3.  0. 

3.  0. 

3.  0. 

3.  0. 

3.  0. 

-.2022E-0!  .  1829E 
>.  0. 

J.  0. 

. 1091E-0I  1789E 


.  1532E-01  0. 


0.  0. 

0.  0. 

0.  0. 

0. 

0.  0. 

0.  0. 

0. 

. 7998E-02  0. 
0.  0. 

0.  0. 

-.5920E-02  0. 
0.  0. 


0. 

0. 

0. 

0. 

0.  0. 

0.  0. 

0.  0. 

.5178E-03  .  1159E 
0.  0. 

0.  0. 

-.290 5E-02  7996E 

0.  0. 


,NFTU,n - IN  FORMAT!  W  1X.E10.9)) 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

1228E-0I  — . 135 IE— 01 

0.  0. 

0.  0. 

.2783E-01  .2I68E-01 

-01  0.  0. 

0.  0. 

0.  0. 

-01  0.  0. 

0.  0. 

0.  0. 

0.  0. 

-.5126E-01  1282E-01 

0.  0. 

0.  0. 

.5016E-01  . 1198E-02 

-02  0.  0. 

0.  0. 

0.  0. 

-02  0.  0. 

0.  0. 

0.  0. 

0.  0. 

-.3112E-01  .1560E-01 

0.  0. 

0.  0. 

. 2022E-OI  -. 1829E-01 
-01  0.  0. 

0.  0. 

0.  0. 

-01  0.  0. 

0.  0. 

0.  0. 

0.  0. 

-.3880E-02  . 1532E-01 
0.  0. 

0.  0. 

-.5178E-03  1 159E-01 

-01  0.  0. 

0.  0. 

0.  0. 

-02  0.  0. 

0.  0. 

0.  0. 

0.  0. 

. 2629E-02  .5920E-02 
0.  0. 


0. 

0. 

0. 

0. 

0. 

0. 

.2853E-02 

0. 


.2972E-02  .92*!7E-02  0. 


-.2972E-02  -.9297E-02 


.37251-02  0. 

0.  0. 


-.2853E-02  -.3725E-02  0. 


FRECNH.FCENT  ZINCH#MFI.A0FP,NH3lNFTlt 
.H00t*09  t200CE*02  !200C4*01  0.  0. 

IN  FOWNTI  3M./,5<  IX,EI0.*l).2t  U.I10)./ .95X1 
<  RWK  I ), I  =1  #N^T& .  1  >  IN  F0WWt<9UX,E10.9>> 

0.  .1006c*01  .2000E+01  .3000001 

.9000E+01  .5000C+01  ,600C€+01  .700(001 

.800CE+01  .  9000C«'0l  .1000002  .110(002 

.  120CE+02  .130(002  .190(002  .150(002 

THE  HEM  FIELD  COMPONENT  IS  EFHI 
IN  FOWINTt 29X.810./.25X) 

CCMNVEEHC  IWW;jC0t)'lH0«=l,m3,l  >.JCOL=l,NFTU,  1 )  — 

0.  0.  0.  0.  0. 

0.  0.  .3730E-01  .2999E-01  0. 

-.3730E-01  -.2999E-01  0.  0.  0. 

0.  0.  0.  0.  0. 

0.  0.  0.  0.  0. 

.  3728E-01  .289CE-01  0.  0.  -.3728E- 

0.  0.  0.  0.  0. 

0.  0.  0.  0.  0. 

0.  0.  0.  0.  . 3710E- 

0.  0.  -.3710E-01  -.2528E-01  0. 

0.  0.  0.  0.  0. 

0.  0.  0.  0.  0. 

0.  0.  .3697E-01  .2039E-01  0. 


IN  F0aP»T(MlX,El0.'O) 

0. 

0. 

0. 

0. 

0. 


-.3728E-03  -.2880E-01 


.  3710E-01 


-.3&97E-01  -.2039E-01  0. 


•2039E-01  0. 


.  3505E-01 


0. 

0. 

. 1938E-0I 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

0. 

0. 

0. 

-.3263E-01 

0. 

0. 

.29I9E-01 


-.3505E-01 


-.7993E-02  0. 


. 3263E-01 


-,2919E-01  -.  1817E-02  0. 


.2993E-0!  -.3921E-02  0. 

0.  0.  0. 

0.  0.  0. 

0.  0.  0. 

0.  0.  -.2018E-0J 

0.  0.  0. 

0.  0.  0. 

0.  0.  . J526E-0I 

J528C-01  . I005E-0I  0. 

0.  0.  0. 

0.  0.  0. 

. I059E-01  -.I139E-01  0. 


.  I817E-02  0. 


0. 

0. 

0. 

-.2993E-01 

0. 

0. 


1938E-01 

0. 

0. 

•  7993E-02 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

.3921E-02 

0. 

0. 


,  7918E-02  0. 


.  J52EE-01  -.1005E-0! 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

-.E390E-02  . 1 169E-01  0. 


.2016E-01  -.7918E-02 

).  0. 

>.  0. 

>.  0. 

).  0. 

>.  0. 

>.  0. 

).  0. 

-.1059E-0!  .II39E-0I 

>.  0. 

).  0. 


.839(002  1169E-01 


Data  Set  II 


.11706-01  -.87276-0*  0.  0.  0. 

0.  0.  0.  0. 

FftEOHH.PCCNT,  ZINCH.AOFl  ,80FP J  W3,  W>TU 
.110«+01  .20006*02  !t000*+01  .30006+02  0. 

IN  FOnmTt  38X,/.5(  IX,610.1).2<  1X.I 10>,/,15X) 

in  fcwtfhmlLtio.ii) 

0.  .10066+01  .20006 +01  . 30066+01 

.10006+01  .50006 ♦Cl  .60006+01  .70006+01 

.80006+01  .10006+01  .10006+02  .11006+02 
.12006+02  .13006+02  .11006+02  .15006+02 

THE  N6M  FIELD  COMPONENT  IS  EPHI 
IN  F0W»7I  29X,A10,/,25X> 

(<U8V6FN(IB0N.JC0L)ltn01»l=l.W3.1).JC0L=l.M»m.l>  -- 


IN  FOBHftK  6(  IX. 610.1) ) 


.16636-01  -. 

.’awr-w4: 

.  1 775E-01  . 

.21176-02  . 
.70056-01 


.3722E-03 
. 13126-07 
-.80566-02  - 
-.2111E-01  - 
.21366-05  - 
.52136-05 
.2657E-01 
. 1279E-02 
.1151E-08  - 
.56056-03  - 
I.  0 

-.56056-03 
-.56756-06 
-.81206-02 
-.21116-01  - 
-.19206-01  - 
. 2971E-01 
.21226-01 
.61526-02  - 
.83816-06  - 
-.18616-03  - 
>.  0 
.18616-03 
-.16106-06 
-.16816-02 
-.11856-01 
-.31566-01 
.18816-01  - 
.76126-02  - 
.12336-03  - 
-.10326-06  - 
-.26176-03  - 
).  0 
.26176-03 
.33666-06  - 
.8H06-03 
-.18526-02 
.56106-05 
-.61356-05  - 


5757E-01 


56156-01 

20106-02 

11226-07 


.16636-01  .5757E-01  0. 

0.  0.  0. 

°ill66£-03  -li>85E-W-! 

0.  0. 

-.11666-03  .37656-01  -. 


I511E-03 

3156E-06 

32716-02 

13506-01 

1665E-01 

2561E-01 

32196-01 

11816-02 

3663E-06 

7052E-03 

70526-03 

15086-06 

21116-02 

,11136-01 

,32636-01 

25266-01 

.33116-02 

37096-02 

,39386-06 

65HE-03 

I65HE-03 
.1832E-06 
.13056-02 
.5121E-02 
.69356-05 
.  1810E-01 
.67376-02 
.31816-02 
.  7015E-06 
.  1261E-01 

1 12616-01 
,  3255E-06 
,  1531E-02 
.60456-02 
,  1156E-01 
.2631E-05 


-.58116-02 

-.2065E-01 

-.12726-05 

.2936E-05 

.21116-01 

.8056E-O2 

-.1312E-07 

.6518E-03 

0. 

-.65186-03 

1967E-06 

-.9271E-02 

-.2691E-01 

-.10226-01 

.19206-01 

.21916-01 

.81206-02 

.56756-06 

.6639E-01 

0. 

-.66396-01 
-.69796-06 
-.3837E-02 
-.16596-01 
-.35506-01 
.31566-01 
.  1 185E-01 
.16896-02 
.16106-06 
-.35106-03 
0. 

.35106-03 
.59616-06 
.6H6E-03 
-.12126-02 
-.  39116-05 
-.56106-05 
.  1852E-02 
-.89106-03 
-.33666-06 
-.59216-01 


-.33316-02  . 

-.5155E-01  . 

-.71686-05  . 

. 1665E-01  . 

.13516-01  0. 
.32716-02  -. 
-.31566-06  -. 
-.52726-03  -. 
0. 

.52726-03  -. 
.21966-06  . 

.26-176-01  . 

-.2169E-01  . 

-.31856-01  . 

.32636-01  . 

. 1112E-01  0. 
-.2111E-02  -. 
-.1508E-06  -. 
-.7962E  03  -. 
0, 

.  7962E-03 
.80286-06  . 
.H59E-02  . 
.2209E-02  . 

-.  1016E-01  . 

-.69356-05 
-.5121E-02  C. 
-.13056-02  . 

-.9832E-06  . 

-.  1321E-03 
0. 

.18216-03 
.  1 333E-06  . 

.21256-02  , 

.67106-02  -. 
.18906-01  -, 
-.11566-01  -. 
-.6015E-02  0, 
153  IE— 02 
.32556-06  - 
, 3135E-09 


29976-02 

17756-0! 

20986-06 

Trm=ur 

20656-01 

58116-02 

1579E-08 

58796-03 

5879E-C3 

1151E-08 

9279E-02 

2657E-01 

52936-05 

10226-01 

2691E-01 

9271E-02 

1967E-06 

31306-03 

3930E-03 
8381E-06 
•  6152E-02 
2122E-01 
.2971E-01 
,35506-01 
,  1659E-01 
,  3837E-02 
,69796-06 
.33786-03 

!  33786-03 
,10326-06 
.12336-03 
.76126-02 
.  1881E-01 
.  3991E-05 
.12126-02 
.69166-03 
.59616-06 
.11316-03 

\  1131E-03 
.10306-06 
.  81 15E-03 


0. 

9:il  BF3T 
-.2090E-02 
-.5615E-01 
-.11586-05 

.5155E-01 
.3331E-02 
-.  12176-06 
-.3262E-03 
0. 

.32626-03 

.36636-06 

-.19816-02 

-.32996-01 

-.25696-01 

.3185E-01 

.2169E-01 

-.26176-01 

-.2196E-06 

-.8092E-03 

0. 

.80926-03 

.39386-06 

.37096-02 

-.33196-02 

-.25266-01 

.10166-01 

-.2209E-O2 

-.91596-02 

-.80286-06 

-.1229E-03 

0. 

.12296-03 

.7015F-06 

.3181E-02 

.67376-02 

.18106-01 

-.18906-01 

-.67106-02 

-.2125E-02 

-.13336-06 

.5U1E-01 

0. 

-.51116-01 
-.38966-06 
.99186  -03 


JCOL 


JCOL 


-o5- 


.27*16-03  -.51516-02  0.  0. 

-.81156-03  -.***-03  .5*2*6-©*  -.3*35E-0*  . 
.10306-04  ,38*66-04  -.36106-06  -.11**6-06 
-.37346-0*  -.*7216-05  .73026-03  .  7*856-03  -. 

0.  0.  .7*826-03  .*33*6-02  -. 

.373* -0*  .*7216-05  .28506-05  -.21*06-05  . 
-.287*6-04  .  205*6-04  .1*026-05  .17506-05  -. 
.74816-03  .58*46-03  -.14516-02  -.35**6-02  0. 
.1*516-02  .  35**6-02  -.74816-03  -.58946-03  . 
-.1*026-05  -.17586-05  .  287*6-06  -.205*6-04 

THE  NEAR  FIELD  COMPONENT  IS  ERHO 
IN  FOAMATl  2*3,810  /  25X1 

U  »JRV6EMt  IRON' JCOC. )’ IR0«=1 , MM3, 1  ).JC0L=1  ,NPTU, 
0.  0/  ’  0.  6.  0' 
0  .  0.  -.5757E-01  .16436-01  0 

-.57576-01  .16636-01  0.  0.  0 

0.  0.  0.  0. 

-.11126-05  .2*076-06  .3791E-0*  .11736-03 
-.56116-01  .16836-01  -.1336E-02  -.32256-02  - 

.21736-02  -.2951E-02  .37*16-0*  .1173E-03  - 

-.108*6-07  .  20*06-0*  -.1032E-06  .*0216-07  - 

.1*6*6-03  .37876-03  .3861E-02  -.55126-02  - 
-.2*676-02  -.56786-02  -.51*16-01  .17326-01 

, l*4*E~03  .3787E-03  -.60166-05  .185*6-05  - 

-.15346-06  .29056-06  -.10846-0*  .**6*6-05 

.*5*86-02  -.72946-02  -.*5286-01  .1785E-01  - 

-.*5286-01  .17856-01  .*5*86-02  -.72*46-02 

-.10846-0*  .**6*6-05  -.15366-04  .2*056-06 

-.12106-0*  .*5116-05  . 37*5E-03  .67306-03 

-.36736-01  .180*6-01  -.3*736-02  -.71816-02  - 

.*07*6-02  -.80*06-02  .37*5E-03  .67306-03  - 
.2*726-06  .733*6-06  .11606-05  .8585E-0*  - 

.3*026-03  .62056-03  .282*6-02  -.7772E-02  - 
-.**28E-02  -.63**6-02  -.272*6-01  .17636-01 

.3*026-03  .62056-03  -.128*6-0*  .*56*6-06 
.19656-05  .23556-06  -.1672E-0*  -.33136-05 

.135*6-02  -.678*6-02  -.18326-01  .16226-01  - 

-.1832E-01  .16226-01  .13596-02  -.678*6-02 

-.16726-0*  -.33136-05  .1965E-05  .2355E-06 

-.21586-0*  -. 1567E-05  .3*126-03  .3*226-03 
-.1089E-01  .13936-01  -.**056-02  -.29*36-02  - 

.83126-0*  -.5*2*6-02  .3*126-03  .3*226-03  - 
.19516-05  -.83956-04  .10526-05  -.15816-05  - 

.32086-03  .2597E-03  -.79836-03  -.3979E-02  - 
-.373*6-02  -.122*6-02  -.5*676-02  .11136-01  - 
.3208E-03  .2597E-03  -.2209E-O*  .5*386-05 

-.1009E-06  -1*756-05  -.15*86-0*  .123*6-0* 

-.122*6-02  -.2627E-02  -.20106-02  .82646-02  - 
-.20106-02  .82666-02  -.12' *6 -02  -.26276-02 
-.15986-0*  .123*6-0*  -.100*6-04  -.  l*75E-05  - 

-.63*06-05  .  1*0*6 -0*  .210*6-03  .  25296-0*  - 

-.17036-03  .  57*86-02  -.16386-02  .  58*16-03  - 
-.12*06-02  -.15306-02  .210*6-03  .25296-0*  - 
-.76526-06  -.70086-06  -.67916-06  .1107E-06 

.12576-03  -.2**7E-0*  -.**076-03  -.78926-03 
-.79826-03  .  4*0*6-03  .56UE-03  .  38536-02  - 
.12576-03  -.2**7E-0*  .13**6-05  .9964E-05  - 

-.162CE-06  .***8E-06  .38156-05  .37**6-05 

-.67166-03  -.*0716-03  .72206-03  .26*06-02  - 

.72206-03  .26*06-02  -.67166-03  -.*0716-03 
.38156-05  .  37**6-05  -.16206-06  .***86-06 

.19716-05  -.*1*86-06  .30786-0*  -.30156-05  - 


-.27*16-03 

.6*356-05 

-.28506-05 

-.7982E-03 

-.73026-03 

.36106-06 

-.553*6-0* 

0. 

.553*6-0* 


.51586-02 

.26316-05 

.2190E-05 

-.*33*6-02 

-.7*856-03 

.11**6-06 

-.1*056-0* 

0. 

.1*056-0* 


,1) - IN  FORMAT!  6<  1X.E10.*)) 

'  0. 

0. 

i.  0. 

•.1089E-07  .20906-09 
.21736-02  -.2*5  IE-02 
-.56116-01  . 1683E-01 

-.11126-05  .  2*076-06 

-.6014E-05  . 1859E-05 
-.51916-01  . 1732E-01 

.38616-02  -.55126-02 
-.10326-06  .*021E-07 

.2837E-03  .60036-03 
-.3328E-02  -.70066-02 
.28376-03  .  60036-03 

.2*726-06  .  733*6-06 
.*07*6-02  -.80*06-02 
-.34736-01  .180*6-01 

-.12106-0*  ,*511E-05 
-.128*6-0*  .*56*E-06 

-.2729.-01  .17636-01 

2829E-02  -.77726-02 
.11606-05  .  8585E-06 

.365*6-03  .51*36-03 
-.*6l5E-02  -.*7*56-02 
.365*6-03  .51*36-03 
.1*516-05  -.83*56-06 
.83126-0*  -.5*2*6-02 
*.  1089E-O1  .13936-01 

-.2158E-0*  -.1567E-05 
-.220*6-0*  .5*386-05 
-.5*4  7E-02  .  U13E-01 
-.7*836-03  -.3979E-02 
.10526-05  -.15816-05 
.2812E-03  . 1282E-03 

-.2717E-02  .*7936-05 
.28126-03  .12826-03 

-.74526-06  -.7008E-06 
-.12*06-02  -.15306-02 
-.17036-03  .57*86-02 
-.63*06-05  .  l*0*E-0* 

.  13**E-05  .99666-05 

.5611E-03  .  38536-02 

-.9S07E-03  -.78926-03 
-.6791E-06  .1107E-06 

.5*286-0*  -.25066-0* 
-.33196-03  .**036-03 

.5*286-0*  -.25066-0* 

.  2607E-06  .  26*26  -06 

-.*3826-03  --2773E-03 


.73WE-0  3  .1154C-02  -.1449E-0J  .2552E-03  . 
-.*3#2E-03  -.2773E-03  . 30?«C-0*»  -.30151-05  . 
.2407E-O4  .2492E-04  .3053E-04  -.8804E-07  -.1 
.3Z58E-0*  .1029E-O*  -.3M3E-03  -.2544E-03  .1 
-.1301E-O3  .U12E-03  .8085E-03  .1538E-02 
.3258E-0*  .1029E-0*  -.833CE-04  -.9838E-04  . 
,T**3E-07  -.259IE-04  -.19ME-05  .5828E-04 

-.33I8E-03  -.2021E-03  .9271E-03  .U48E-02 
.927IE-03  .U48E-02  -.3318E-03  -.2021E-O3  J 
-.in*E-05  .5828E-04  .7453E-07  -.2591E-04 

the  nem  field  cwtoneni  is  ezee 

IN  FOWWTt  29X,M0,/.25X) 

<  <  M8VEE7N  I ww’ JCOL ) ,IR0W=1  ,MH3, 1 ) .  JC0L=1  ,NPTO, 
0.  0.  0.  0.  0. 

0.  0.  0.  c. 

0.  0.  0.  0.  0. 

0.  0.  0.  0. 

.3592E-04  .3*HE-07  -.6589E-C5  -.1409E-0* 

.479CE-02  -.4392E-02  -.1396E-01  -.1085E-O1  . 

-.m8€-0*  .5238E-03  -.6589E-05  -.1409E-0*  . 

.3715E-08  .3755E-08  .909*E-07  .7753E-07  . 

-.4001E-0*  -.1077E-03  H29E-03  .1957E-02  . 

-.1272E-01  -.8203E-02  .1297E-0I  -.12*8E-01 

-.6001E-0*  -.1077E-03  .9299E-0P  -.5047E-07  . 

.*392£-04  .2751E-04  .1317E-0*  -.2N20E-O5 

-.3543E-03  .3916E-02  .1431E-01  -.1792E-01 
.  1631E-01  1792E-01  -.3543E-03  .3916E-02 

.1317E-0*  -.2920E-05  .5392E-06  .2751E-04  . 
.2IH5E-0*  -.8981E-05  -.9253E-03  -.3882E-03 
.  177*E-OI  -.22*51-01  -.880CE-02  -.9509E-03  . 
-.3163E-03  .5W7E-02  -H253E-03  -.3882E-03  . 

.%**E-04  .3*52E-04  .1208E-O5  .5005E-07  . 

-.6250E-O3  -H022E-03  .2459E-03  .7210E-O2  . 
-.4367E-02  .2295E-02  .1639E-01  -.2585E-01  . 

-.42*€C-03  -.*022E-03  .2339E-0*  -.1487E-0*  . 

.  8523E-06  -.3572E-06  .196GC-0*  -.215IE-0*  -. 

. 1331E-02  .  7708E-OZ  .1235E-01  -.2782E-01 

. I235E-01  -.2782E-01  .133IE-02  .7708E-02  -. 

.1940E-0*  -.2151E-0*  .8523E-04  -.3572E-04  . 

.  1Y52E-0*  -.2205E-0*  -.7303E-03  -.1798E-03  . 

.4423E-02  -.2805E-01  -.  1W5E-02  .5691E-02  . 

.  259*C-02  .  7355E-02  -.7303E-03  -.1798E-03  . 

.1339E-04  -H040E-06  -.3858E-04  .1353E-07  . 

-.4477E-03  -.5251E-0*  .3742E-02  .6351E-02  . 

. 1338E-02  .5*77E-02  .3052E-03  -.26*HE-01  . 

-.4677E-03  --5251E-CH  .  1051E-0*  -.2093E-O* 

-.3450E-O4  .5*28E-04  .6849E-05  -.2012E-O* 

,*435£-02  .*932E-02  -.5578C-02  -.2317E-01  . 

-.5570E-O2  -.2317E-01  .5435E-02  .N932E-02 

.6849E-0S  -.2012E-O*  -.3450C-04  .5*28E-04  . 

.2397E-05  -.19I6E-0*  -H805E-03  .1397E-03  . 

-.1027E-01  1843E-01  .3977E-02  .2810E-02  -. 

.5095E-02  .3319E-02  -.5805E-03  .1397E-03  . 
.5209E-07  .7*H3E-04  .5378E-04  .5058E-04  -. 

-.37*6€-03  . 2028E-03  .5103E-02  .1722E-02  -. 

,*UlE-02  . 1275E-02  -.1333E-01  -.1339E-OI  . 

-.37*«-03  . 2028E-03  -.2M1E-05  -.1656E-0*  . 

H841E-04  .8302E-O7  -.625 *€-05  -.U97E-0*  -. 

.97J4E-02  .322SE-03  -.H48E-01  -.8079E-O2  . 

-.H48E-01  -.8079E-O2  .5716E-02  .3229E-03 

-.625*E-05  1197E-0*  H841E-04  .8302E-O7  , 

711  IE-05  -.6652E-05  -.1750E-03  .2361E-03  . 


.7345E-03 
.  1971E-05 
-.8330E-04 
.80  85E-O3 
•.3513E-03 
•  3053E-04 
.*0*1E~0* 
■.9345E-0* 
HOME-0* 


.1954E-02 
-.9H8E-06 
.9838E-06 
.1530E-O2 
-.2*64E-03 
-.8806E-07 
.  3502E-05 
.  I529E-03 
,3*02E-05 


—  in  formakmjxfio  *■> 
0. 

,I**OE-Ol  -.1186E-0I 
0. 

.  3715E-08  .3755E-08 

.  1558E-0*  .5238E-03 

•  6790E-O2  -.6392E-02 
.  3592E-06  . 3* l*E -07 

H259E-05  -.5067E-O7 
.  1297E-01  -.12*8E-01 
.1929E-03  .1957E-02 

.909*E-07  .7753E-07 

.  2061E-03  -.2636E-03 
.  1095E-O1  -H660E-02 
.2061E-O3  -.2636E-03 
.  9&**E-04  ,3*52E-06 

.3J63E-03  .5897E-02 

.  1779E-01  -,22*^-01 
.2H5E-0*  -.8981E-05 
.2339E-0*  -.16C7E-0* 

.  163*E-01  -.2585E-01 
.  265*£-03  . 7210E-02 

.  '208E-O5  .50WE-07 

.7286E-03  3126E-03 

.3721E-02  H582E-02 

.7286E-03  -.3126E-03 
.  1339E-06  -H060E-O6 
.2595E-02  . 7355E-02 

.6623E-02  -.2805E-91 
.1552E-0*  -.2205E-0* 

.  1051E-0*  -.2093E-0* 

.  3052E-03  -.26«WE~01 
.  3762E-02  .  6351E-02 
.3858E-06  . 1353E-07 

.5791E-03  .5370E-O* 

.3078E-02  H373E-02 
.5791E-03  .5370E-O* 

.5209E-07  .7553E-06 

.5095E-02  .3319E-02 

.I027E-01  1863E-01 

.2397E-05  -.  1916E-0* 

.2691E-05  -. 1656E-0* 

.  1333E-01  -.  1339E-01 
.5103E-02  . 1722E-02 
H378E-06  .5W8E-06 
.  2676E-03  .2352E-03 
.  3762E-C2  .5651E-0* 

.2676E-03  .2352E-03 

.2362E-06  -.1750E-06 
H068E-02  -.7675E-03 


OO  J  © 


-.19596-01  -.3Z92E-02  .3219E-02 
.S06HE-O2  -.76756-03  -.17906-03 
.2362E-06  -.17906-06  -.38376-07 
-.10596-03  .21806-03  .33006-02 
.Z635E-02  -.19286-02  -.13326-01 
-.10596-03  .21806-03  -.56296-05 
-.11326-06  .62536-07  -.35816-05 
.25256-02  -.20386-02  U*fO€-01 

-.11906-01  .90326-02  .2525E-02 

-.3581E-05  -.19096-05  -.1132E-06 
THE  NEAR  FIELD  COMPONENT  IS  HPHI 
IN  FORMAT!  29X,A10,/,25X> 

C ( WAVEEW I ROM, JCOC ) , I R0W=1  ,MM3, 1 
0.  '  0.  '  ' 
0.  -.59826-01 

•5982E-01  .2175E-01  0. 

0.  0. 

. 10736-05  -.18266-06  .1737E-03 
-.5320E-01  .2275E-01  .16656-02 

. 9989E-03  -.96126-02  .1737E-03 
.98336-07  . 19536-08  .5895E-06 
.5532E-03  . 1157E-03  .12606-02 
.29216-02  -.39736-02  -.98526-01 
.5532E-03  .1157E-03  .9591E-05 
.  1890E-O5  -.51906-06  .  23506-05 
.5981E-03  10I1E-01  -.91386-01 

-.9138E-01  .27396-01  .5981E-03 

. 235CE-05  -.93066-05  .1890E-O5 

-.12296-09  -.2995E-05  .11UE-02 

-.3281E-01  .2870E-OI  .33206-02 

-.6969E-03  -.10586-01  .HUE-02 

. 307TE-05  -.18386-05  .2907E-05 
.12096-02  .29966-03  -.22096-02 
.29966-02  -.73906-02  -.2387E-01 
.12096-02  .  29966-03  -.31396-09 

. 1197E-05  -.96986-05  -.39336-09 
-.35796-02  -.8972E-02  -.I539E-01 
-.15396-01  .28036-01  -.35796-02 

-.3933E-09  .  33306-09  .U97E-05 

-.28176-09  .5353E-C9  .10936-02 

-.7661E-02  .  26206-01  -.5007E-03 
-.9592E-02  -.6603E-02  .10936-02 
-.1172E-05  -.9221E-05  -.2687E-05 
.  7879E-03  -.3820E-03  -.99666-02 
-. 1993E-02  -.60996-02  -.1078E-02 
.  7879E-03  -.38206-03  -.39206-05 
-.27206-05  -.2792E-06  .1991E-09 

-.98806-02  -.2868E-02  .92636-02 

.92636-02  .2001E-01  -.98806-02 
.  1991E-09  .  97336-09  -.27206-05 

.30926-09  .2523E-09  .23286-03 

.82586-02  .1597E-01  -.3837E-02 
-.99506-02  1975E-02  .2328E-03 

-.1609E-05  .10726-05  -.2681E-06 

.  7791E-09  -.3116E-03  -.3892E-02 
-.90516-02  -.18596-02  .10896-01 

.  7791E-09  -.31166-03  .  27006-09 

.52206-06  .63306-06  .15966-09 

-.33686-02  .19396-03  .12236-01 
.12236-01  .7369E-02  -.33686-02 

.1596E-09  -.61936-05  .52206-06 

.95296-05  -.6770E-O5  .3310E-09 


-.8069E-03 

.2361E-03 

-.13376-06 

-.1535E-02 

.8195E-03 

-.3216E-05 

-.1909E-05 

.9032E-02 

-.2038E-02 

.62536-07 


-.19596-01 
-.711 1E—05 
-.56296-05 
-.13326-01 
.33006-02 
-.3837E-07 
-.61636-09 
•  2022E-02 
-.61636-09 


-.32926-02 
-.6852E-05 
-.32166-05 
•  8195E-03 
-.15356-02 
--1337E-06 
. 1991E-03 
-.18996-02 
.1991E-03 


).JCOL=l,NPTU,l) 
0.  0. 


IN  FORMAT!  6(  1X,E10. 9)) 
0. 

.2175E-CI  0.  0. 

0.  0.  0. 

0.  .9833E-07  .19536-08 

.2961E-09  .99896-03  -.96126-02 
-.2080E-02  -.53206-01  .2275E-01 

.2961E-09  .10736-05  -.18266-06 

-.9728E-07  .9591E-05  -.16716-05 

-.81506-02  -.98526-01  .2506E-01 
.2506E-O1  . 1260E-02  -.81506-02 

-.  1671E-05  .58*56-06  -.9728E-07 

-.9306E-O5  .8986E-03  .  22996-03 
.27396-01  .3501E-02  -.55366-02 

-.10116-01  .89866-03  .22996-03 

-.5190E-06  . 3077E-O5  -.1838E-05 
.29916-03  -.69696-03  -.10586-01 
-.66806-02  -.3281E-01  .2870E-01 

.2991E-03  -.12296-09  -.2995E-05 
-.3617E-05  -.31396-09  .10736-09 
-.9912E-02  -.23876-01  .28906-01 

.28906-01  -.2209E-02  -.99126-02 
.  1073E-09  .2907E-05  --3617E-05 
.33306-09  .11866-02  .67626-09 

.28036-01  .10796-02  -.7966E-02 

-.89726-02  .11866-02  .6762E-09 
-.96986-05  -.11726-05  -.9221E-05 
-.17836-03  -.95926-02  -.66036-02 
-.7037E-02  -.7itiE-02  .26206-01 
-.17836-03  -.2817E-09  .53536-09 
-.23996-05  -.39206-05  .5929E-09 
-.96386-02  -.  1078E-02  .  23996-01 
.2399E-01  -.9966E-02  -.96386-02 
.5929E-09  -.26876-05  -.2399F-05 
.97336-09  .  98896-03  -.96686-03 
. 2001E-01  -. 31506-02  -.9786E-02 
-.28686-02  .98896-03  -.96686-03 
-.27926-06  -.16096-05  .10726-05 
-.9259E-03  -.99506-02  -.1975E-02 
-. 3302E-02  .82586-02  .15976-01 
-.92596-03  .  3092E-O9  .2523E-09 

. 1252E-05  .27006-09  .97136-05 

-.98826-03  .10896-01  .1165E-01 

.  1166E-01  -.38926-02  -.98826-03 
.9713E-05  -. 2681E-06  .1252E-05 

-.61936-05  .  29976-09  -.1985E-03 

.73696-02  -.38896-02  -.58596-03 
.19396-03  -  2997E-09  -.19856-03 

.63306-06  .  55966-06  -.9389E-07 

-. 1337E-03  -.2929E-02  .72086-03 


-38- 


O I |  t I  |  O I  | I  ©  o© 


.12966 -01  .33566-02  -.39756-02  .99996-03 
-.29296-02  .7 2086-03  .33106-09  -.13376-03 
.55966-0*  -.93896-07  .19956-06  -.93196-06  - 
*SJ|“09  ^.12l?E-03  -.25086-02  . 1 196*32 
-.29056-02  .1291E-02  .11786-01  -.2t9C6-03  - 
“.12196-03  69806 -06  -.21936-05 

-.29356-0*  -.32116-0*  -.3091E-06  .18826-05 
-.20356-02  16296-02  .10356-01  -.32316-02  - 

.10356-01  -.32316-02  -.2035E-02  .16296-02 
-.30916-06  18826-05  -.29356-06  -.32UE-06 

the  neo*  field  component  is  hr  ho 

IN  FO«mT(29X,ftlO./,25X) 

<  ( W9VE6W  IRON, JCOt  l^IROWs^ma,  1 ) . JC0L=1,NPTU; 
.  0.  0.  u.  0; 

-  _  0.  -.21756-01  -.59826-01  0 

.21756-01  .59826-01  0.  0.  0 

0.  0.  0. 

.21936-06  .  97936-06  -.3026E-09  .16286-03 
.22626-01  -.51866-01  0.  0. 

.99136-02  -.93906-03  .  30266-09  -.16286-03  - 
.97986-09  -.96816-07  .86276-07  .5131E-06 

.12266-03  .92696-03  .6807E-02  .9835E-03  - 

0.  .22756-01  .93736-01  - 

.12266-03  -.92696-03  -.22796-05  -.29226-05  - 
.70256-06  .13106-05  .53026-05  -.62706-05  - 

.66176-02  .9509E-03  -.18706-01  -.33026-01  0 
.18706-01  .33026-01  -.6697E-02  -.95096-03 

.53026-05  .  62706-05  -.70256-06  -.13106-05 

.58966-07  -.25086-09  -.29196-03  .  9157E-03 

.10996-01  -.23986-01  0.  0. 

.97806-02  -.19196-03  .29196-03  -.9157E-03  - 

.19576-05  -.1229E-05  .27166-05  -.9530E-06  - 

.18686-03  .2711E-03  .23706-02  .95096-03  - 

0.  .31696-02  .17626-01  - 

.18686-03  -.27116-03  .18136-09  .3396E-09  - 

.  1679E-05  -.2632E-05  -.35396-09  -.17086-09 
.27616-03  .13936-02  .18836-02  -.15386-01  0 

-.18836-02  .15386-01  -.2761E-03  -.13936-02  - 

.35396-09  .17086-09  -.16796-05  .26326-05  - 

-.31226-09  .16196-09  .12996-03  -.25686-03  - 

.39036-02  -.1515E-01  0.  0. 

.93256-03  -.28236-02  -.1Z99E-03  .25686-03 
.89196-06  .33596-05  -.31916-05  -.18086-05  - 
.98086-05  -.55596-03  -.90296-03  .9367E-02 
0.  0.  -.2195E-02  .15286-01 

-.98086-05  .  55596-03  .19066-05  -.3922E-09 

-.35016-05  .99276-06  .352;r’-09  .32966-09  - 

.35986-03  .53936-02  -.53996-03  -.1981E-01  0 

.53996-03  . 1981E-01  -.35986-03  -.5393E-02 

-.35216-09  -.32966-09  .  35016-05  -.99276-06  - 

.53396-09  .17826-05  -.6277E-03  -.9891E-03 

-.38086-02  -.13396-01  0.  0. 

-.23026-02  -.52906-02  .6277E-03  .9891E-03  - 

.I698E-05  -.28996-05  .  75106-06  .  28126-05 

-.7561E-03  -.12126-03  .90296-02  .3910E-02  - 

0.  0.  .68196-02  .10726-01  - 

.75616-03  .12126-03  -.91886-09  .31286-09  - 

.21036-05  .11066-05  .11706-09  -.99286-09  - 

.97206-02  .17726-02  -.88196-02  -.7157E-02  0 

.88196-02  .71576-02  -.97206-02  -.17726-02 

-.11706-09  .  99286-09  -.21036-05  -.11066-05 

-.19976-09  -.3285E-09  -.33156-03  .39766-03 


.12966  -01 
.95296-05 
-.69886-06 
.11786-01 
-.25  08E-O2 
. 1995E-0b 
.92666-09 
-.22376-02 
.92666-09 


.33566-02 
-.67706-05 
-.21936-05 
-.2190E-03 
.11966-02 
-.93196-06 
-.  1 37  IE— 03 
.18166-02 
- .  1 37  IE— 03 


,1)  -  IN^FORMATCM  1X,E10.9>) 

'!  o! 

o. 

.97986-09  .9681E-07 

•  9913E-02  .93906-03 

•2262E-01  .51866-01 

.21936-06  -.97936-06 
.22796-05  .29226-05 
.2275E-01  -.93736-01 
.6807E-02  -.98356-03 
-.  8627E-07  -.5131E-06 
.29326-03  .50326-03 

0. 

.29326-03  -.50326-03 
.19576-05  .12296-05 

.97806-02  . 1919E-03 

.10996-01  .23986-01 

.58966-07  .  25086-09 

. 1813E-09  -.3399E-09 
.31696-02  1762E-01 

.23706-02  -.95096-03 
.27166-05  .95306-06 

. 1139E-09  .59036-09 

0. 

.11396-09  -.59036-09 
.89996-06  -.33596-05 
.93256-03  .28236-02 
.39036-02  . 1515E-01 

.  3122E-09  -.16196-09 
. 1996E-05  .39226-09 

.  2195E-02  -.15286-01 
. 90296-03  -.93676-02 
. 3191E-05  . 1808E-05 
.30396-03  -.66326-03 
0. 

.30396-03  .6632E-03 
.16986-05  .  28996-05 

.23026-02  .52906-02 

.38086-02  .13396-01 

.5339E-09  -.17826-05 
.91886-09  -.31286-09 
.68196-02  -.10726  -01 
.90296-02  -.39106-02 
.  75  1  06  -06  -.2812E-05 
.62156-03  .23296-03 

0. 

•  6215E-03  -.2329E-03 

.17076-05  -.71916-06 
.92006-02  2231E-03 


«.  ’ll  -  *  T. 1*^*  «J' 


%  -•  v*  r*  vr' 


-.92631-02  -.32751-02  0.  0. 

-.92001-02  .22311-03  .3315E-03  -.39761-03  . 

-.17071-05  .7191E-04  .30261-06  -.13921-05 

•*.70171*09  .39961-03  .29601-02  -.13761-02  -. 

0.  0.  .81891-02  .5768E-05  -. 

.70171-09  -.39961-03  .  23061-09  .10921-09  -. 
-.82131-06  -.78081-06  -.19201-09  .59761-05  . 

.17661-02  -.16131-02  -.63031-02  .19981-02  0. 

.63031-02  -.19981-02  -.17661-02  .16131-02  -. 
.19201-09  -.59761-05  .82131-06  .78081-06 

THE  TEA*  FIELD  COMPONENT  IS  HZEE 
IN  FORMAT!  29X,A10./,25X) 

t  (  WAVEEHC  I  ROW,  JCO£. ) ,  I  R0W=1  ,MM3, 1 )  JC0L=1 .  NPTU , 
0.  0.  0.  6.  0. 
0.  0.  0.  0.  0. 
0.  0.  0.  0.  0. 
0.  0.  0.  0. 

.78921-07  -.55901-06  .28731-06  -.69091-09  -. 

.98361-02  . 10971-01  0  .  0. 

.16001-02  . 1699E-03  -.28731-06  .69091-09  *. 

-.9132E-08  . 1352E-07  .59651-07  -.3281-1-06  . 

.  1259E-09  -.9059E-03  -.53581-02  -.1555E-09  . 
0.  0.  -.22211-01  -.1739E-01  . 

-.12591-09  .90591-03  -.91921-07  .55191-05 

-.35061-07  -.  1572E-05  -.2966E-05  -.10691-09  . 

-.91131-02  . 1311E-02  .3999E-01  .18031-01  0. 

-.39991-01  -.18031-01  .91131-02  -.13111-02 

. 2966E-05  .10691-09  .35061-07  .15721-05  -. 

-.7668E-05  .12781-06  . 1667E-03  -. 12221-02  -. 

.9157E-01  . 1238E-01  0.  0. 

.11291-01  -.36931-02  -.16671-03  .12221-02  . 

.11191-05  .39091-05  -.3533E-05  -.91091-05  -. 

.19181-03  -.1998E-02  -.U22E-01  .63791-02  . 

0.  0.  -.91931-01  -.39121-02  . 

-.19181-03  .19981-02  .2677E-06  -.29931-09  . 

-.59251-05  -.2299E-05  .2959E-09  .5999E-09  . 

-.92951-02  .86931-02  .35561-01  -.51991-02  0. 

-.35561-01  .51991-02  .9295E-02  -.86931-02  -. 

-.29591-09  -.59991-09  .59251-05  .2299E-05  -. 

.69811-09  .9919E-09  -.38101-03  -.19921-02 

.26931-01  -.11181-01  0.  0. 

.60071-02  -.98891-02  .38101-03  .19921-02  -. 

.6299E-05  -.13001-05  -.38121-05  .93561-05  . 

-.8096E-03  -.11361-02  -.21771-02  .9695E-02  . 

0.  0.  1801E-01  .19091-01  . 

.80961-03  .11361-02  -.9268E-09  -.1051E-O9  . 

-.36291-07  .  9875E-05  .79UE-09  -.37711-09  -. 

.11671-02  . 7857E-02  .10231-01  -.19221-01  0. 

-.10231-01  .19221-01  -.11671-02  -.78571-02  . 

-.79111-09  . 3771E-09  .36291-07  -.98751-05  . 

.35891-09  -.69521-09  -.88891-03  .1207E-09  . 

,92*9E-02  -.12981-01  0.  0. 

-.31231-02  -.5068E-02  .8889E-03  -.1207E-09  -. 

-.27131-05  -.27851-05  .30021-05  -.17201-06  -. 

-.5121E-03  . 3965E-03  -33761-02  -22951-02  . 

0.  0.  -.39731-03  .96761-02  -. 

.51211-03  -.39651-03  .99701-05  .56121-09  -. 

.12851-05  -.19991-05  -.32881-09  -.29291-09  -. 

.29131-02  .98981-03  -.19121-02  -.6819E-02  0. 

.19121-02  .68191-02  -.29131-02  -.98981-03  . 

.32881-09  .  29291-09  1285E-05  .19991-05  -. 

-.26991-09  .61531-05  .9665E-09  .16631-03  . 
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.  32851-09 
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-.57681-05 
.13761-02 
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0. 

0. 
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98361-02  -.  1097E-01 
78921-07  .  '590E-06 

91H2E-07  -.55191-05 
2221E-01  . 1739E-01 

,53581-02  . 1555E-09 

5965E-07  . 3285E-06 

70931-09  -.86631-03 
0. 

,70931-09  .86631-03 

,11191-05  -.39091-05 
,11291-01  .36931-02 

,91571-01  -.12381-01 
7668E-05  -.12781-06 
2677E-06  .29931-09 

9193E-01  .39121-02 

1122E-01  -.6379E-02 
35331-05  .9109E-05 

,10701-09  -.1567E-02 
0. 

107CE-O9  . 1567E-02 

,62991-05  .13001-05 

.60071-02  . 98B9E-02 

,26931-01  .11181-01 

•  6981E-09  -.99191-09 
,92681-09  .10511-09 

,  1801E-01  -.19091-01 
,21771-02  -.96951-02 
.38121-05  -.9356E-05 
,10161-02  -.5605E-03 
0. 

,10161-02  .  5605E-03 

,27131-05  .  27851-05 

,31231-02  .50681-02 

,92991-02  .12981-01 

.35891-09  .69521-09 

,99701-05  -.56121-09 
,39731-03  -.96761-02 
.33761-02  -.22951-02 
.30021-05  .17201-06 

,  1169E-03  . 3603E-03 

0. 

,11691-03  -.36031-03 
,  7263E-06  -.17321-05 
,12091-02  -.565 6E-09 
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-.1815E-02  -.9762E-02  0. 
-,1209E-02  .5656E-09 
.T263E-06  .  1732E-05 

.8593E-09  -.3222E-09  . 

0.  0. 

-.8593E-09  .3222E-09  . 
-.9989E-0*  .9977E-0<>  . 

.6221E-03  .5995E-03 
. 2172E-02 


•H762E-02  0.  0.  .1815E-02 
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.6153E-05 
. 1731E-09 
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.9136E-05  .  9729E-05  -.3830E-04  -.9682E-09 


.2172E-02  -.3302E-02  0.  0. 
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THE  NEAR  FIELD  COMPONENT  IS  EPHI 
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.30&5E-O1  -.2680E-01  0. 

0.  0. 

.3759E-07  . 1052E-06  .1820E-09 

- 3065E-01  -2589E-01  0. 

. 1519E-02  .2220E-03  -.1820E-O9 

-.1298E-09  -.3098E-10  .1236E-08 

.6282E-09  -.5790E-09  -.2866E-02 
0.  0.  -.3056E-01 

-.6282E-09  .5790E-09  -.2831E-06 

.9913E-09  -.2751E-08  .8913E-06 

-.3906E-02  1561E-09  .3015E-01 

3015E-01  -.1913E-01  .3906E-02 

-.8913E-06  -.2091E-O5  -.9913E-09 
.1991E-05  . 3819E-05  .1369E-03 

.2919E-01  .1A02E-01  0. 

.95I6E-02  -.5216E-03  1369E-03 

.I081E-07  .1028E-07  -.9912E-07 

.1299E-03  -.2520E-03  -.9692E-02 
0.  0.  -.2739E-0I 

-.1299E-03  .2520E-03  -.3669E-05 

-.9191E-07  .2695E-09  .5772E-05 

-.A308E-02  .2009E-02  .2968E-01 

-.2968E-01  -.3180E-O2  .9308E-02 

--5772E-05  -.5776E-05  .9191E-07 

. 7789E-05  .9823E-05  .2560E-09 

.2132E-OI  W9E-02  0. 

.3610E-02  -.2633E-02  -.2560E-09 
.1265E-06  -.99I3E-07  -.1201E-06 
-.9079E-09  -.2839E-03  -.2697E-02 
0.  0.  1799E-01 

.9079E-09  .2839E-03  -.8922E-05 

--7115E-07  . I323E-06  .8619E-05 

1729E-02  .3120E-02  .1399E-01 

I399E-01  . 7557E-02  .1729E-02 

-.8619E-05  .2960E-O6  .7U5E-07 

.6929E-05  -.2556E-05  -.1237E-03 
. 9619E-02  899^-02  0. 

.695*€-03  -.2967E-02  .1237E-03 
.9523E-08  -.1197E-06  .2678E-07 


.2680E-01  0.  0. 

I.  0.  0. 

>.  .129BE-09  . 3098E-10 

•.I501E-09  -.1519E-02  -.2220E-O3 
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.2329E-01  .2866E-02  .2635E-03 
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.1891E-03  1991E-05  -.3819E-05 
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.  1298E-02  .2739E-01  .8502E-02 

•8502E-02  .9692E-02  -.1298E-02 
.5288E-05  .9912E-07  .I509E-07 

•  5776E-05  .8833E-09  -.2991E-03 

. 3I80E-O2  0.  0. 

.2009E-02  -.8833E-09  .2991E-03 
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,  9952E-07  . 8922E-05  .  255SE-05 

.3019E-02  .1799E-01  -.5069E-O2 

.5069E-02  .2697E-02  -.3019E-02 
.255  9E-05  . 1201E-06  -.9952E-07 

.2960E-06  -.9302E-O9  -.2357E-03 
.7557E-02  0.  0. 

.3I20E-02  . 9382E-09  .2357E-03 

.  1323E-06  95  23E-08  .1197E-06 

.  1761E-03  -.8959E-03  .2967E-02 

-.9619E-02  . 8999E-02 

.  1761E-03  -.6929E-05  .2556E-05 

. 7J27E-07  .9561E-05  -.3660E-05 


IX.  MICROFICHE  OF  THE  LISTING  OF  THE  FORTRAN  SOURCE  OF  PROGRAM 
PCOEF  AND  LISTINGS  OF  THE  AZIMUTHAL  COEFFICIENTS _ 

The  microfiche  is  48X.  There  are  a  total  of  thirty-four  fiches. 
Thirty- three  are  listings  of  the  azimuthal  coefficients  that  were 
written  on  the  nine  track  tapes.  The  other  fiche  is  the  Fortran 
Source  listing  of  the  Program  PCOEF.  Program  PCOEF' s  source  was 
also  included  on  one  of  the  nine  track  tapes.  It  can  be  used  to 
read  the  azimuthal  coefficients  from  the  tapes. 

The  azimuthal  coefficients  are  for  a  mine  buried  in  one  of  three 
soil  types  (silt  loam  either  at  5%,  10%  or  20%  water  content)  where 
eleven  frequencies  are  used  for  the  incident  field  (400,  500,  600, 

700,  800,  900,  1000,  1100,  1200,  1300  and  1400  Mega  Hertz).  Each 
fiche  contains  the  listing  of  the  azimuthal  coefficients  for  one  of 
the  thirty-three  combinations  of  soil  type  and  frequency. 

Each  fiche  comes  in  its  own  package.  All  packages  are  labeled 
as  well  as  each  fiche. 

The  package,  for  the  fiche  containing  the  Fortran  Source  listing 
of  the  Program  PCOEF,  is  labeled  "SOURCE".  The  fiche  is  labeled 
"SOURCE  OF  PCOEF". 

The  label  on  the  package  for  the  fiche  containing  the  azimuthal 
coefficients  consists  of  the  letter  "C"  followed  by  six  digi  The 
first  four  digits  are  the  frequency  in  Mega  Hertz.  The  last  two  digits 
are  the  water  content  in  percent  for  the  soil.  Thus,  the  label 
"C040005"  is  for  400  Mega  Hertz  silt  loam  having  5%  water  content. 

While  the  label  "Cl 30020"  is  for  1300  Mega  Hertz  silt  loam  having 
20%  water  content.  The  fiche  containing  the  azimuthal  coefficients 
is  labeled  with  the  frequency  and  the  water  content  in  percent.  For 
example  the  fiche  for  the  azimuthal  coefficients  corresponding  to  an 
incident  field  of  400  Mega  Hertz  for  silt  loam  having  5  percent  water 
content  is  labeled  "400  MHZ  05  PERCENT";  while  the  one  for  1300  Mega 
Hertz  20%  water  content  is  labeled  "1300  MHZ  20  PERCENT". 

For  a  detailed  description  of  the  azimuthal  coefficients  see 
Chapter  VII  "Description  of  the  files  containing  Azimuthal  Coefficients 
(sent  via  the  nine  track  tapes)".  For  a  description  of  the  Program  PCOEF 
see  Chapter  VI  "Description  of  the  Program  PCOEF  (sent  via  nine  track 


X.  CONCLUSION  AND  RECOMMENDATION- 


The  Unimoment  method  was  used  previously  to  calculate  the 
scattered  fields  by  a  buried  land  mine  18].  There  was  some  doubt  on 
the  credibility  of  the  results  because  scattered  fields  were  found 
to  be  stronger  at  lower  frequencies,  contrary  to  the  experiences  of 
MERADCOM  scientists.  By  using  a  more  realistic  model  of  the  ground 
parameters,  we  have  indeed  discovered  that  the  scattered  fields  of 
the  dielectric  mine  buried  in  silt  loam  having  a  water  content  of  20% 
are  stronger  at  1000  MHz  than  at  either  600  or  1400  MHz  which  is 
consistent  with  the  experiences  of  MERADCOM  scientists.  This  result 
suggests  that  the  operating  frequency  of  the  existing  mine  detectors 
have  been  appropriately  chosen. 

The  massive  amount  of  calculated  results  obtained  in  this 
work  deserves  attention,  since  it  contains  information  that  may  be  useful 
for  the  improvement  of  mine  detectors.  We  recommend:  1)  that  a  program 
to  study  these  results  be  initiated  to  utilize  it  for  engineering 
purposes  and  2)  that  an  experimental  program  to  verify  the  data  be 
initiated. 

The  present  results  are  limited  to  mines  which  are  buried  at 
least  1.6"  deep.  In  warfare,  antipersonnel  land  mines  are  buried 
almost  flush.  The  Unimoment  method  may  be  modified  to  compute  flush 
buried  land  mines.  It  is  expected  that  a  shallow  buried  land  mine  will 
yield  greater  scattering  than  a  deeply  buried  mine,  thus  enhancing  the 
probability  of  detection  from  shallowly  buried  mines.  A  third  recommenda¬ 
tion  is  that  the  present  computer  program  should  be  generalized  to  compute 
electromagnetic  scattering  from  flush  buried  antipersonnel  land  mines. 
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A  primary  field  is  defined  as  the  sum  of  the  incident  field 
and  the  reflected  field  from  the  interface,  without  the  scatterer. 
Well  known  formulas  of  primary  fields  are  available  (for  example 
pages  299-301  of  reference  Ill]).  For  the  sake  of  completeness 
we  list  the  fields  of  both  polarizations  in  MKS  units  as  follows: 

(i)  The  case  of  magnetic  field  perpendicular  to  the  plane  of 
incidence,  (0°  polarization  angle): 


Referring  to  Figure  1-1,  the  incident  field  is 
..Inc  _  -1  Ax  Az 


_i  JBvx  jB  z 
—  e  x  e 
% 


where  3X  =  kQ  sine,  3Z  =  kQ  cose,  k0  =  “  * 

X  -  c  =  2.99793  *  108  meters/sec 
o  f 

The  primary  fields  for  z  >  0  are: 
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<y  Vector  points  out  of  page 

S)  Vector  points  into  page 

The  y  axis  points  into  the  page  and  so  the  coordinate 
system  is  right  handed. 


Figure  1-1 


Field  orientations  of  0°  and  90° 
Polarization  incident  waves 
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(iii) 
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The  primary  field  for  the  case  of  6=0°,  polarization 
angle  =  0°,  for  1400  MHz,  at  Z  =  1  in.  when  the  soil 
is  silt  loam  having  a  water  content  of  20%  as  used  in 
the  sample  data,  may  be  found  as  follows: 


eg  =  15.6  -  jl.869  (From  Table  2  on  page  12) 


29.34  k  Z  =  .7453  (in  MKS  units) 
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=  95.043  ej*0596  ohms 
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0° 


03„x  r  j3  z  -OB  z' 
e  +  pe 


-1 


3Z  =  kQ  cos  30°  =  25.410  [meter  J 


3zZ  =  .64541 
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(A17) 


B .  Example  of  Computing  the  Scattered  Fields  for  0°  Incidence 


Using  Set  I  of  the  file  in  Chapter  VIII  (polarization  angle 
0°)  EPHI  component  at  Z  =  1  in.  can  be  computed  as  follows: 

(i)  To  compute  EPHI  for  4>  =  0°,  (positive  x-axis) 
at  p  =  0,  add  the  11  complex  numbers  in  block 
labeled  JCOL  =  1.  The  summation  of  the  block 
labeled  JCOL  =  2,  gives  EPHI  for  <{>  =  0° , 
p  =  1  cm,  etc.  For  this  particular  incident 
plane  wave  polarization,  EPHI  =  0  on  the 
x-axis  as  to  be  expected. 

(ii)  To  compute  the  same  field  for  <}>  f  0°,  say 
<f>  =  30°,  we  use  the  formula. 


M 


-M 


Since  this  is  the  case  of  0°  incidence,  only 
am  with  m  =  ±1  are  non-zero, 


For  p  =  0, 


^  e^30 

9  -1 


*  ♦  a.,^30' 


*+r 


=  (.02202  +  j. 06676)  e"j30°  -  (.02202  +  j. 06676)  eJ’30° 


=  -2 j(. 02202  +  j. 06676)  sin30° 


=  .06676  -  j.02202,  for  6  =  0°,  <j>  =  30°,  p  =  0°  (Bl) 
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C.  Example  of  Computing  the  Scattered  Fields  for  30°  Incidence 


R 


Using  Set  II  of  the  file  in  Chapter  VIII,  the  scattered  EPHI 
component  can  be  computed  as  follows: 

(i)  To  compute  EPHI  for  <j>  =  0°,  (positive  x-axis) 
at  p  =  0,  add  the  11  complex  numbers  in  block 
labeled  JCOL  =  1.  The  summation  of  the  block 
labeled  JCOL  =  2  gives  EPHI  for  4>  =  0°, 
p  =  1  cm  etc.  For  this  particular  incident 
plane  wave  polarization,  EPHI  =  0  on  the 
x-axis  as  to  be  expected. 

(ii)  To  compute  the  same  field  for  <j>  =  0°,say 
<j>  =  30c,  we  use  the  formula 

5  ■  £  4  e« 

-M 


for  p  =  0,  (only  am  with  m  =  ±1  are  non-zero) 


’4> 


=  a,  g-j308  +  a  .  eJ*30' 
-1  +1 


=  (.01663  +  j. 05757)  e"j30°  -  (.01663  +  j. 05757)  ej30° 


=  -2 j(. 01 663  +  j. 05757)  sin30° 


=  .05757  -  j. 01663  for  6  =  30°,  Q  =  30°,  p  =  0 


L 
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(Cl) 


p  =  1  cm  (use  data  in  block  JCOL  =  2  in  set  II) 


=  (-.7005  x  10"9  +  j.1122  x  10"7)  e"J’5  x  30° 
+  (.2098  x  10"6  +  0.1158  x  10"5)  e~J*4  x  30° 

+  (.1166  x  10"3  -  j.3785  x  10“4)  e“J'3  x  30° 

+  (-.2997  x  10"2  -  j.2090  x  10"2)  e"j2  x  30° 

+  (.1775  x  10_1  +  j.5615  x  10"1)  e"j30° 

+  0 

-  (.1775  x  10" 1  +  0-5615  x  10"1)  ej30° 

-  (-.2997  x  10"2  -  j.2090  x  10“2)  eJ’2  x  30° 

-  (.1166  x  10~3  -  j.3785  x  10"4)  e°’3  x  30° 

-  (.2098  x  10“6  +  j.1158  x  10“5)  eJ'4  x  30° 

-  (-.7005  x  10“9  +  j.1122  x  10“7)  eJ'5  x  30° 


Since  the  coefficients  for  m  4  are  negligible  in  comparison 
to  the  lower  order  ones,  we  shall  omit  them  in  the  following 
summation,  which  gives: 

E*  =  -[  2j(.  1166  x  10"3  -  vi. 3785  x  ICf4)  sin90° 

+  2j(-  2997  x  10"2  -  j  .2090  x  10"2)  sin60° 

+  2j ( .  1 775  x  10"1  +  j.5615  x  10"1)  sin30°] 

=  -£(.0000757  +  j. 0002332)  +  (.003619  -  j. 0051 90) 

+  (-.05615  +  j.01775)J  =  .05245  -  j. 01279 

for  0  =  30°,  <j)  =  30°,  p  =  1  cm  (C2) 


D.  Example  of  Computing  Total  Fields  for  0°  Incidence 


Using  the  file  of  Chapter  VIII,  the  total  EPHI  component 

at  z  =  1  in.  can  be  computed  as  follows: 

✓ 

(i)  To  compute  total  EPHI  for  <J>  =  0°,  (positive  x-raxis) 
at  p  =  0,  simply  add  the  primary  field  to  the 
scattered  field  computed  in  (B).  Using  formula  (A6) 
at  =  0°.,  we  find  EPHI  primary  to  be  zero.  From 
Section  B(i)  the  scattered  field  for  this  component  is 
zero.  Thus,  as  expected,  the  PHI  component  of  the 
total  electric  field  is  zero  along  the  x-axis. 

(ii)  To  compute  the  same  field  for  $  f  0°,  say  <j>  =  30°, 
we  add  the  primary  field  to  the  scattered  field 
already  computed  in  (B).  Using  the  formula  (A2), 

(A6)  and  (A17), 

E?  =  (-.3089  -  j  1.0974)  sin<J> 

<p 

=  -.1545  -  j . 5487  for  6  =  0°,  <J>  =  30° 

The  primary  field  is  independent  of  p  because 
8  =  0°. 

Using  Bl,  we  find  the  total  field  is 
E^  =  E*  +  Ej  =  (.06676  -  j. 02202)  +  (-.1545  -  j.5487) 
so 

E.  =  -.08774  -  j . 5702  =  .5769  e"0*1’723 


In  the  time  domain 


E<j>(t)  =  Re^5769e"j1'723  ejtotj  =  .5769  cos  (<ot  -  1.723} 
for  0  =  0°,  <j>  =  30°,  p  =  0° 

Using  (Dl)  and  (B2)  we  find 

E<*>  =  Zl  +  =  (*0^498  -  j. 02340)  +  (-.1545  -  j .5487) 

=  -  .08952  -  j. 57210  =  .5791e“J’ 1,7260 
for  0  =  0°,  <J>  =  30°,  p  =  1  cm 


(D3) 


(04) 


In  the  time  domain 

ty*)  =  Re^5791e"j1,7260  eJcotj=  .5791  cos  (wt  -  1.7260) 
for  0  =  0°,  4>  =  30°,  p  =  1  cm 


(05) 
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E.  Example  of  Computing  Total'  Fields  for  30°  Incidence 


The  total  EPHI  component  at  Z  =  lin.  at  30°  incidence  can 
be  computed  as  follows: 

(i)  To  compute  total  EPHI  for  <p  =  0°,  (positive  axis)  at 
p  =  0,  simply  add  the  primary  field  to  the  scattered 
field  computed  in  (C).  Using  formula  (A6)  at  4>  =  0°> 
we  find  EPHI  primary  to  be  zero.  From  section  C(i),  the 
scattered  field  for  this  component  is  zero.  Thus,  the 
PHI  component  of  the  total  electric  field  is  zero  along 
the  x-axis. 

(ii)  To  compute  the  same  field  for  <J>  f  0°.  say  4>  =  30°, 
we  add  the  primary  to  the  scattered  field  already 
computed  in  (C).  Using  the  formulas  (A2),  (A6),  and 
(A18), 

Ej  =  -.8828  eJ* 11,1991  +  14,671x1  sin<J> 

<P 

and  since  x  =  0  at  p  =  0 

Ej  =  -(.3206  +  j.8226)  sin  30°  =  -.1603  -  j.4113  =  -.4414e3’ 

for  0  =  30°,  <p  =  30°,  p  =  0 
for  P  f  0,  we  use  x  =  p  cos<j> 

eJ  =  -.8828  1,1991  +  J’14-671(p  cos  30°)^  sin  30° 

♦ 

=  -.8828  1,1991  +  *1271'  sin  30° 

=  -.8828  e° 1,3262  sin  30°  =  -.1069  -  j . 4283 


.1991 

(El) 


for  0  =  30°,  <f>  =  30°,  p  =  1  cm 


(E2) 


Using  (Cl)  and  (El)  we  find  that  the  total  field 

=  E*  +  Ej  =  (.05757  -  j. 01663)  -  (.1603  +  J .4113) 

=  -.1027  -  j . 4279  =  .4401e'° 1,806 
for  6  =  30°,  =  30°,  p  =  0 

In  the  time  domain 

E^t)  =  Re|*4401e~^ 1,806  eJw^=  .4401  cos(wt  -  1.806) 

for  6  =  30°,  (t>  =  30°,  p  =  0 

Using  (C2)  and  (E2)  we  calculate  for  the  total  field 

E^  =  E*  +  E^  =  (.05245  -  j. 01279)  -  (.1069  +  j.4283) 

=  -.05445-  J .4411  =  .4444e"J*1,69/! 
for  6  =  30°,  0  =  30°,  p  =  1  cm 

In  the  time  domain 

E.(t)  =  Rf.4444e->i  1-694  gO^t'l  =  >4444  c6s(wt  _  1>694) 


z 


Appendix  II 


Source  of  PCOEF 


0 


B 


■ 
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ooooooooooo  ooo  OO  O  O  ooooo  O  o  O  O  OOOOOOOO  OO  OOOOOOOOOO 


PROGRAM  FC0EF(  IN^UT,  rAPElGr.INPUTtTAPEl9» 
1  OUTPUT, TAPE16=0UTPUT) 


THE  PURPCSE  OP  THIS  PROGRAM  IS  TO  PRINT  THE  AZIMUTHAL  COEFFICIENTS 
THAT  ARE  GOING  TO  BE  DELIVERED  TO  THE  ARMY. 


THE  MAIN  PART  OF  T^IS  °RCGRAM  WAS  LAST  CHANG-  0  ON 
-  i*  JANUARY  1932  - 


COMPLEX  A ZIC (15,16,6) 


DIMENSION  RHO  < 16 ) ,RH0C<16) ,AC0MP(6> ,WCCMF(o ) ,  RC ONP ( 6) , M C 15 ) 


NOTE  ONLY  UP  TO  THE  FIRST  15  AZIMUTHAL  COEFFICIENTS  FOR  EACH 
COMPONENT  AT  EACH  RHO  Z  POINT  PAIR  CAN  SE  USED,  ALSO  NOTE  THAT 
THAT  TAPE,  ITAPfcl*  MUST  HAVE  EXACTLY  16  RHO  VALUES  FOR  EACH  7. 
FURTHER  MORE  THEY  MUST  3E  APPROXIMATELY  AS  LISTED  IN  THE  ARRAY, 
RHOC# 

DATA  IR71/»ITA.°EI/l8/,MDIM/i5/,MAXFTS/i6/*  NPTSC/16/  » IFLAGE/1/ 

DATA  RBOC/O.tl.,2.,3.,4.,5.,6.,7.,3.,9. ,1Q.,11«,  12»,13*,14»,15»/ 

THE  ALLOWED  NAMES  OF  THE  FIELD  COMPONENTS  WILL  NCW  BE  SET, 

DATA  ACOM  P/13H  EPHI  *10H  EPHO  , 1 CH  EZEE 
1  10H  FT  A0*HPHI » 10  H  ET  AO*HRHO»  10H  E7AG*HZEF/ 

THE  NAMES  OF  THF  FIELD  COMPONENTS  WHICH  HILL  BE  USED  BY 
SUBROUTIWE  WAZIC  IN  WRITING  THE  RESULTS  IN  A  HUMAN  REACA9LE 
FORM  WILL  NOW  St  SET . 

DATA  WCONP/1QH  EPHI  ,10H  ERHO  ,1GH  EZEE  , 

1  10H  HPHI  ,1  OH  HRHO  , 1  OH  HZEE  / 

OAT A  IUN ITR/10 / , IUNIT E/16/, IUN IT0/1 6/ 

DATA  ISET/0/ 


1  CONTINUE 

THE  DESIRED  COEFFICIENTS  HAVE  9CZN  OBTAINED  FLAG  WILL  9F  UNSET. 
IFOES= ISET+1 

SUBROUTINE  READIN  READS  THE  INPUT  FILE.  THE  INPUT  FILE 
HAS  THE  INFORMATION  CONCERNING  THE  FREQUENCY  IN  MEGA  HERTZ, 

THE  WATER  CONTENT  OF  THE  SCIL  IN  »£R  CENT,  THE  ANGLE  OF  TNCinp.NCE 
IN  DEGREES  OF  THE  INCIDENT  PLANE  WAVE,  AS  WELL  AS  THE  ANGLE  OF 
POLARIZATION  IN  DEGREES  CF  THE  INCIDENT  PLANE  WAVE.  ALSO,  The; 
DISTANCE  IN  INCHES  ABOVE  THE  GROUND  WHERE  THE  NEAP.  FIELDS  APE 
EVALUATED  IS  RE  AO »  SUBROUTINE  READIN  READS  THE  NUMBER  CF  OIFFERFNT 
COMPONENTS  OF  T He.  NEAP.  FIELD  THAT  ARE  TC  BE  PRINTED  IN'  AODITION  TO 
THE  NAMES  OF  THE  COM3ONENTS  IN  THE  ORDEP  THEY  ARE  TO  BE  PRINTED. 


CALL  REACIN (IUNIT5* IUNITR, 

1  FMHZW ♦ FCENTW*AOFIW,AOFPW,Z INCHW , NCGMP .RCCMP) 

SUBROUTINE  COHCHK  CHECKS  TO  3c  SURE  THE  COMPONENT  NAMES  ARE 
ALLOWED. 

CALL  COMCHK  (IUN  IT  E,A  COMP,  NCOMP,  ROOM  P> 

SUB ROUTT N  E  RW1W1Z  MAKES  THE  FIRST  READ  TO  T APE  I *  THE  TAPE  WHERE 
THE  AZIMLTHAl  COEFFICIENTS  ARE  STOREO.  SU3 ROUT IME  RWlWiZ  WILL 
READ  IN  THE  AZIMUTHAL  COEFFICIENTS  OF  THE  NEAR  FIELD!  THAT 
CORRESPOND  TO  ONE  INCIDENT  PLANE  WAVE. 

CALL  RWiMZ(IR»ITAaEI*IUNITE» 

1  MAX  FT  S»  FREOMH,  PCENT , Z I NCH , A CFI , AOFf  , 

2  MUSED  ,  NPTSJ  » RKG,  MDIM, AZIC) 

SUBROUTINE  TCHECK  MAKES  SURE  THAT  THE  USER  HAS  SPECIFIED  A 
FREQUENCY  AND  A  PERCENTAGE  WHICH  IS  IN  THE  FILE. 

CALL  TCH£ CK (IUNIT E»FMHZWfFR£QMH,PCE NT W, PCENT) 

2  CONTINUE 

SUBROUTINE  IaZCHK  OE I ERMINES  IF  THE  LAST  REAO  OF  TAPE,  TAPEI, 
READ  THE  DESIRED  COEFFICIENTS. 

CALL  IPZCHK(AOFIW,AOFI, AOFPW,ACFP, ZINCHW, ZINCH, 

1  ISET.IFDES) 

IF  IFOES  IS  SET  THEN  THE  LAST  READ  DID  03TAIN  THE  DESIRED 
COEFFICIENTS. 

IF  ( IFDcS  ot.Q.  ISET)  GO  TO  4 

3  CONTINUE 

THE  LAST  READ  WAS  NOT  FOR  THE  DESIRED  C OEFFICEMS. 

SO  ONE  MCRE  READ  WILL  3E  ATTEMPTED. 

CALL  RWiWlZdR,  17 A°EI,IUNIT£, 

1  MAXFTS,  FREQMH,  FCENT, ZI NCH , A OFI , AC FF , 

2  MUSED  ,  NPTS’JtRHO,  MDIM,AZIC> 

GO  TO  2 

4  CONTINUE 

THE  LAST  READ  WAS  FOR  THE  DESIRED  COEFFICIENTS. 

SU3ROUTINE  RHOCHK  MAKES  SURE  ThAT  THE  RHO  VALUES  ARE 
AT  LEAST  APPROXIMATELY  AS  EXPECTED. 

CALL  RHGCHKCUNITE,NPTSU  ,RHO,  NFTSC,RHOC) 

SUBROUTINE  LOADM  LOADS  THE  ARRAY  M  WITH  THE  AZIMUTHAL  CROERS 
OF  THE  COEFFICIENTS. 

CALL  LCACMdUNITE,  MU3E0,  M) 

EACH  LOOP  TO  5  CONTINUE  WILL  PRINT  THE  AZIMUTHAL  COEFFICIENTS 
OF  THE  SERIES  WHICH  REPRESENTS  ONE  OF  THE  COMPONENTS  OF  T'HE 
SCATTERED  FIELD. 

DO  5  1  =  1 , NCOMP, 1 
P.C 0  M P£  =  R  C  OH?  ( I ) 


SU8K0UTISE  CNUME  WILU  RETURN  THE  COMPONENT  INDEX. 

CALL  CNU  h  S ( IUNITE  » RC  3MPE , ACOMP. I  COMP) 

SUBRGUTIhE  WAZIC  WILL  PRINT  THE  COEFFICIENTS. 

IFLAG=IFLAGE 
IFdCOhF  .GT.  3)  IFLA!G=IFLAG  +  1 
CALL  WAZICdUNITO  tMDIM,  RHO*  IFLAG. MUSED, 

1  M,  AOFI.AOFP.  Z  INCH,  RCOMPE»WCCMP(  ICOMP),  AZICd.it  ICOMP)) 

5  CONTINUE 

6  CONTINUE 


OOO  oooo  OO  OOOOOOOOOCJO 


SUBROUTINE  CNUMS(  I'JNITE,  RCOMPE,  ACOPP.ICOHP) 


THIS  SUB  SOUTINE  SETS  ICO.MP  EQUAL  TO  THE  INDEX  NUN8ER  OF 
THE  ELEMENT  OF  ARRAY,  ACOMP,  THAT  IS  EQUAL  TC  RCOMPE. 

IF  RCOMPE  IS  NOT  EQUAL  TC  ANY  OF  THE  ARRAY  ELEMENTS  THEN 
THE  USER  IS  SO  INFORMED  AND  THE  PROGRAM  IS  STOPPED. 


THIS  SUBROUTINE  MAS  LAST  CHANGED  ON  5  NOVEMBER  1961. 


DIMENSION  ACOMP  (6 ) 


1  CONTINUE 

FOR  THE  SAKE  OF  EFFICIENCY  SOME  LOCAL  VARIABLES  WILL  SE  SET  EQUAL 
TO  SOME  CF  THE  PARAMETERS. 

IUNIT=  IUMTE 
RCLOC=RCCMPE 
DO  2  1=1, 6,1 
INDEX= I 

IF  ( RCL  CC  «EQ.  ACOMP  ( I) )  GO  TO  T 

2  CONTINUE 

RCOMPE  DCES  NOT  MATCH  ANY  OF  THE  ALLOWED  VALUES. 

WRITE  ( I'JMT,  1C  C  2)  (ACCNP(  V  ,1=1,6, 1)  ,RCLOC 
103  2  FORMAT (1  hi, 

11H  , 25HTF IS  IS  SUBROUTINE  CNUM5., /, 

21H  , 29HA  f  ERROR  HAS  BEEN  DISCO VERED . , / , 

31H  , 43HRC  GMPE  DOES  NOT  EQUAL  ANY  OF  THE  ALLOWED  VALUES*, 2 (1H  ,/), 
A1H  ,  11  HA COiMP (1 )  =  , A 1 Q,  1  3H  ACOMP(2>  =  »Al0tl3H  AC0MP(3>  =  ,Aio 
5  , 

61H  ,  11 H  AC  OMp  (4  )  =  ,AiG,l3H  ACOMP(5>  =  ,A10,13H  ACOMP  (6  )  =  ,  A  10 

7  ./, 

81H  , 11H  RCOMPE  =  ,A1G,2(1H  ,/), 

91H  ,A3hTHS  SUBROUTINE  WILL  NEXT  STOP  THE  PROGRAM.) 

STOP 

3  CONTINUE 
ICOMP=INCEX 
RETURN 

END 


OOOC5  ooo  nooo  OO  OOOOOOOOOOOOO 


SUBROUTINE  COMCHK  < IUNITE  ,  ACC MP,  NCOFP,  ROOM3) 


THIS  SUBROUTINE  CHECKS  THE  FIRST  NCOMP  ELEMENTS  OF  THE  ARkAY, 
RCTMP,  COMPARING  THEM  TO  THE  ALLOWABLE  SETINGS.  IF  FACE  CF  T’-^E 
FIRST  NC CMP  ELEMENTS  EOUAL  ONE  OF  THF  ALLOWABLE  VALUES  THEN  THIS 
SUBROUTINE  TAKES  NO  FUTHER  ACTION.  IF  JUST  ONE  OF  THE  FIRST 
NCOMP  ELEMENTS  OF  RCOMP  DOES  NOT  EQUAL  ANY  Or  THE  ALLOWABLE  VALUE 
THIS  SUBROUTINE  SO  INFORMS  THE  USER  AND  THEN  STOFS  THE  PROGRAM. 


THIS  SUBROUTINE  WAS  LAST  CHANGED  ON  5  NOVEMBER  1961, 


DIME NS  ION  RCOMP (NCOMP) ,AC0MP  (6 ) 


1  CONTINUE 

FOR  THE  SAKE  OF  EFFICIENCY  SOME  LOCAL  VARIABLES  WILL  BE  SET  FCUAL 
TO  SOME  CF  THE  PARAMETERS. 

IUNIT=  IUN  ITE 
NUMB=NCCMP 


EACH  OF  THE  FIRST  NCOMP  ELEMENTS  OF  RCOMF  WILL  NCW  BE  CHECKEC* 

DO  7  1  =  2,  NUMB,  1 
RCOMPI  =  RC  CMP ( I ) 

DO  2  J=l,6,i 
ACOMPJ=ACOMP(J) 

IF  (  RCO  HP  I  cEQ.  ACOMPJ)  GO  TO  f, 

2  CONTINUE 

3  CONTINUE 


THE  ITH  ELEMENT  OF  THE  ARRAY ,  RCO-'P,  DOES  NOT  EOUAL 
ANY  OF  THE  ELEMENTS  OF  THF  ARRAY,  ACOMP. 


WRIT£(IUMT,10Q3)  I,  ( ACO HP (K ) , K=1 , l) 

1003  FORMAT  (1H* 

11H  *  26HTF IS  IS  SUBROUTINE  COMCHK.,/, 

21H  ,29HA)  ERROR  HAS  3EEN  OISCO VEREO . * /, 

31H  ,  19HTFC  ELEMENT  NUMBER  ,Ii,26H  OF  ARRAY,  RCOHF 
A1H  ,  55HEGLAL  ANY  OF  THE  ALLOWEC  VALUES  PASSED  IN 
5  2  ( 1 F  ,/), 

61H  ,  11HAC CMP (1 )  =  , AiO.,  1  3K  ACCMP(2)  =  ,A10,13H 

7  *  / » 

31H  » 11  HA  C  CMP ( 4 )  =  *Ai0*13H  AC0MP(5)  =  *Al0*13H 
L  CONTINUE 


,  DOEb  NO  I  ,/, 
ARRAY,  ACOMP,, 

ACCMP  (3  )  =  ,  A  IS 

ACC  HP  (6  )  =  ,  A  10) 


DO  5  K=i , NUMB, 1 
WRITE(IUNIT,1004)  K,RCOMP(K) 

100+  FORMAT  f 

11H  ,  FhRCOMPUll,  4H)  =  ,A10) 

5  CONTINUE 

WRITE(IUMT,10G5) 

1005  FORMAT  (IF  ,/, 

11H  ,  A3  HTF  IS  SUBROUTINE  WILL  NEXT  STOP  THE  PROGRAM.) 
STOP 

5  CONTINUE 
C 

C  THE  ITH  ELEMENT  CF  THE  ARRAY,  RCOMP*  EQUALS  THE  JTH 
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'  *  f  — 


ELEMENT  CF  THE  ARRAY,  ACOMP. 

7  CONTINUE 

EACH  OF  THE  FIRST  NCOM=>  ELEMENTS  OF  THE  ARRAY,  P.COMP, 
EQJAL  ONE  CF  THE  ALLOWABLE  VALUES. 

RETURN 

END 


n 


r 
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SUBROUT  I h  E  ECFILE  ( IUNir  .  ISET.  IENCF) 


THIS  SU3R OUTINE  CETERMINES  IF  CN  A  PREVIOUS  READ  OF  UNIT 
THE  ENC  CF  THE  FILE  WAS  REA C  FED# 


THIS  SU3RGUTINE  WAS  uAST  CHANGED  ON  24  JANUARY  1981. 


1  CONTINUE 

IFCEOF  (UNIT))  2.3.2 

2  CONTINUE 

PREVIOUS  READ  ENCOUNTERED  ThE  END  OF  THE  FILE. 

IENOF=ISET 

3  CONTINUE 
RETURN 


SUBROUTINE  IPZCHK < AOFIW, AOFIG,  AOFPR.AOFPG, ZINCHW. ZINCHG, 
1  ISET ♦ IFLAG) 


THIS  SUBROUTINE  COMPARES  AOFIW  TO  AOFIG,  AOFPW  T<~  AOFPG,  ANG 
ZINCHW  TC  ZINCHG.  IF  ALL  COMPARISONS  ARE  CLOSE  THEN  IFLAG  IS 
SET  EQUAL  TO  ISET.  OTHERWISE  THE  VALUE  OF  IFLAG  IS  NOT  CHANGED. 


THIS  SUBROUTINE  WAS  LAST  CHANGED  ON  5  NOVEMBER  1991. 


DATA  smL/i.E-3/ 


1  CONTINUE 

AO  SI C=  ABS  (AOFIW-AOFIG) 
IFtABSIO  .GT.  SMALL)  GO  TO  5 

2  CONTINUE 

A8SPD= ABS (AOFPW-AOFPG) 

IF  (ABS  PC  .GT.  SMALL)  GO  TO  5 

3  CONTINUE 

A8SZ  0= ABS  (ZINCHW- ZINCHG ) 

IF ( A BS  ZC  .GT.  SMALL)  GO  TO  5 
-  CONTINUE 

ALL  COMPARISONS  ARE  CLOSE. 

IFL AG= IS  ET 
5  CONTINUE 
RETURN 
END 


L 
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OOOOOOO  OOOO  OO  OOOOOOOOOCJ 


SUBROUTINE  LOADM( IUNTTE, MUSED, M) 


THIS  SUBROUTINE  SETS  THE  ELEMENTS  OF  THE  ARRAY  M  TO  THE  AZIMUTHAL 
ORDERS  OF  THE  COEFFICIENTS  THAT  WILL  BE  PRINTED.  I  Ft Y  WILL  EE  IN 
THE  ORDER  THEY  ARE  TO  3E  PRIN.ED. 


THIS  SUBROUTINE  WAS  LAST  CHANGED  ON  17  NOVEMBER  1991. 


DIMENSION  M( MUSED  ) 


1  CONTINUE 

FOR  THE  SAKE  OF  EFFICIENCY  SUMS  LOCAL  VARIABLES  WILL  BE  SET  EOUAL 
TO  SOME  CF  THE  PARAMETERS. 

IUNIT  =  11  NI7.E 
MTOTAL  =  ML'SED 

MUSED  IS  THE  NUM3ER  AZIMUTHAL  ORDERS  USED.  THIS  NUMBER  SHOULD  BE 
ODD  SINCE  IT  IS  ASSUMED  THAT  THE  SMALLEST  ORDER  USED  IS  THE 
NEGATIVE  OF  THE  LARGEST  ONE  USED  ANO  THAT  THE  ALL  THE  ONES  IN 
THESE  LIMITS  ARE.  USED  AS  WELL.  IT  IS  ALSO  ASSUMED  THAT  THE 
AZIMUTHAL  ORDERS  ARE  IN  ACSENDING  ORDER. 

MP0S=MT0TAL/2 
MCHECK=MTOTAL- <2*MP0S) 

IF  ( MTOT  Al  *GT  •  C  .ANJ.  MCHECK  .EQ.  1)  GO  TO  3 

2  CONTINUE 

WRITE(IUNIT»10C2)  MT  OTA  L 
1032  FORMAT  llH, 

11H  , 25HTH IS  IS  SUBROUTINE  LOADM.,/, 

21H  *  29HA  f  ERROR  HAS  BEEN  DISCOVERED.,/, 

31H  ,39HMl’SED  SHOULD  BE  A  POSITIVE  ODD  INTEGER. ,  2  (1H  ,/), 

41H  ,  9H  MUSED  =  ,  I10.2(1H  ,/>, 

51H  ,  A3 HT M  IS  SUBROUTINE  WILL  NEXT  STOP  THE  PROGRAM.) 

STOP 

3  CONTINUE 
MVALUE=-M  POS 

DO  4  IM=l .MTOTA  L,  1 
M(IM)=M VALUE 
MVALUE  =  M\(AlUE  +  1 

4  CONTINUE 

5  CONTINUE 
RETURN 


n  n  a  n  oo  ooooooooo 


SUBROUTINE  RCNAME  (IUNITR  ,NUM,RCOMP> 


THTS  SUBROUTINE  READS  FROM  UNIT,  IUNITR,  THE  NAMES  CF  THE 
COMPONtNTS  TO  3f  PRINTED. 

THIS  SUBROUTINE  WAS  LAST  CHANGED  ON  13  NOVEMBER  iP8i. 


DIMENSION  RCOMP  (6  ) 


1  CONTINUE 

FOR  THE  SAKE  OF  EFFICIENCY  THE  UNIT  NUMBER  PASSED  AS  A  PARAMETER 
WILL  8E  SET  TO  A  LOCAL  VARIABLE • 

IUNIT=  IUNITR 

READ  (IUNIT,  1001)  (ROOM9  C II  ,I  =  1,N'JM,1) 

1001  FORMAT (6A10) 

2  CONTINUE 
RETURN 
END 
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no  noo  ooooo  ooo  o  o  o  oooo  oo  o  o  ooooooooooooo 


SUBROUTINE  READIN  CU.NITE,  IUNITR, 

1  PRiQIH,  PCF.NT,  A  OF  I,  A  OF  P,  ZINCH,  NCCMP, RCC^P ) 


THIS  SU  3R  CUT IN£  READS  FRCM  UNIT,  IUNITR,  THE  INF  CRHAT I  ON  NEEDED 
TO  S  FECI  f  Y  WHAT  THE  USER  WANTS  PRINTED.  THE  INFORMATION  IS  THE 
FREQUENCY  IN  MEGA  HERTZ,  THE  PERCENTAGE  CF  ThE  WATER  CONTENT  THAT 
IS  IN  THE  SOIL,  THE  HEIGHT  ABOVE  THE  GROUND  IN  INCHES  THAT  THE 
NEAR  FIELDS  ARE  TO  3E  PRINTED,  THE  NUMBER  OF  DIFFERENT  COMPONENTS 
TO  BE  PRINTED,  AND  THE  NAMES  OF  THE  COMPONENTS  TC  EE  PLOTTED, 


THIS  SUBROUTINE  WAS  UAST  CHANGED  ON  13  NOV^'.ER  1981. 


DIMENSION  PCOMP  (5 ) 


DATA  ISET /0/ , IUNS  ET/1/ 


1  CONTINUE 

FOR  THE  SAKE  OF  EFFICIENCY  THE  UNIT  NUMBERS  PASSED  AS  PARAMETERS 
WILL  BE  SET  TO  LOCAL  VARIABLES. 

IERRsIUMTE 
IREAD=  IUMTR 

THE  FOUNT  END  OF  FILE  FLAG  WILL  NOW  3E  UNSET. 

IENDFs IUN  SET 

THE  FIRST  READ  OF  UNIT,  IUNITR,  WILL  NOW  BE  GONE. 

READ  (I READ,  1001)  FMHZ^C,  AI,AP,  Z,RNUM 
1031  FORMAT (IX ,F9.4,5F10.4> 

IF  THE  ENO  OF  FILE  MARK  WAS  READ  BY  THE  LAST  REAT  A  SERIOUS  ERRCR 
HAS  BEEN  DISCOVERED.  SUBROUTINE  EOF ILE  DETERMINES  IF  THE  PREVIOUS 
READ  OF  UNIT,  IUNITR,  ENCOUNTERED  THE  £NC  OF  FILE  HARK • 

CALL  E OF ILE ( IRE AC , ISET, I ENO F ) 

IF  (ISET  .NE.  IENDF)  GO  TO  3 

2  CONTINUE 

THE  ENO  CF  FILL.  MARK  WAS  ENCOUNTERED  DURING  THE  FIRST  READ 
CF  UNIT,  IUNITR,  THE  USER  WILL  PE  SO  INFORMED  AND  THEN  THE 
PROGRAM  Y  ILL  BE  STOP°ED* 

WRITEt IER  R«  103  2  )  I  REA  D 
1002  FORMAT  (1  H» 

11H  , 26 HTp IS  IS  SUBROUTINE  READIN.,/, 

21H  , 29HA  L  ERROR  HAS  BEEN  DISCOVERED. , /» 

31H  ♦  33  HRE  ACHED  THE  END  OF  THE  FILE  ON  UNIT  =  ,13,/, 

41 H  ,  25  HA  F  TER  ONLY  THE  FIRST  READ.,/, 

51H  ,43HTF:S  SUBROUTINE  WILL  NEXT  STO^  THE  PROGRAM . ) 

STOP 

3  CONTINUE 

NUMs  IF  IX  (  RNU  F ) 


c 

c 


NUM  SHCULO  BE  A  POSITIVE  INTcGIR  LESS  T HAM  OR  LQUAL  TO  6 


C 

C 

C 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


I F  ( 0  •  LT  •  NUM  .AND.  NUM  ,LE.  6)  GO  TO  5 
4  CONTINUE 

AN  ERROR  HAS  3EFN  DISCOVERED. 


WRITE  ( IEF  R,  I<j3  4)  RNUM,NUM,FMHZ,  PC  ,A  I*  A  D,  7. 

100  4  FORMAT  (IK, 

11H  t  2&  HTF IS  IS  SUBROUTINE  REAOIN.,/, 

21H  »  29HA  f  ERROR  HAS  SEEN  DISCOVERED../. 

31 H  .  5SHP.COMp  SHOULD  BE  A  POSITIVE  NUMBER  LESS  THAN  CF  EOJAL  TO 
4  ,2(1H  ,/), 

51H  ♦  9H  RNUM  =  .EiG.L.iiX*  9H  NCOMP  =  .110  »/, 

61H  ,  9HFFFQMH  =  . E10.4.11K.  9H  PCENT  =  .E1C.4,/. 

71H  ,  9H  AOFI  =  .E10.4.11X,  9H  AOFP  =  ,E10.4,/, 

81H  .  9H  ZINCH  =  ,Ei£.4,2(iH  ./). 


=  .E1G.4,  UK, 
=  »  E10 . 4. 1 IX  . 
=  .E10.4.11X, 
=  ♦ EiC .4,2 (iH 


♦110  ./. 
♦E1C .4,/. 
» El  0.4,/, 


29HNCCMP  WILL  BE  SET  EQUAL  TO  1.) 


NUM  =  i 
CONTINUE 

SUBROUTINE  RON  A  ME  WILL  NOW  READ  FROM  UNIT,  IUNITF.  T^e  FIELD 
COMPONENTS  THE  USER  WANTS  PRINTED. 


RCN A  ME ( IREAD »NU M,  F.CCMP) 


IF  ON  THE  LAST  FE  AD  GF  THh  UNIT ,  IJNITR.  THE  ENO  OF  FILf  MARK  WAS 
READ  THEN  A  SERIOUS  MISTAKE  HAS  DCCUREO.  SUBROUTINE  EOFILE 
DETERMINES  IF  THE  °RFVIOUS  READ  OF  UNIT.  IUNITR,  ENCOUNTCREO  THE 
END  OF  FILE  MARK. 


CALL  EOFlLEdREAC.TSET,  EENDF) 

IFUSET  .  NL.  It NCF)  GO  TO  7 
5  CONTINUE 

THE  END  Cf  FILE  MARK  WAS  ENCOUNTERED  DURING  THl  SECOND  READ 
OF  UNIT.  IUNITF.,  THE  USER  WILL  BE  SO  INFORMED  AND  THEN  THE 
PROGRAM  k  ILL  BE  STOPPED. 

WRITE ( IEF  p, 100  6 ) IREA  D 
iOOb  FORMAT  (lUi 

11H  .26HTHIS  IS  SUBROUTINE  REAOIN.,/, 

21H  » 29HA  f  ERROR  HAS  BEEN  DISCOVERED.,/. 

31H  .  38  HRE  ACHED  THE  END  OF  THE  FILE  ON  UNIT  =  ,13,/, 

41H  ,  27HAFTER  ONLY  THE  SECOND  READ.,/, 

51H  » 43  FT  F  IS  SUBROUTINE  WILL  NEXT  STOP  THE  PRCGFAM.) 

STOP 

7  CONTINUE 

DEMONSTRATE  THAT  NOTHING  ELSE  READ  FROM  THE  UNIT,  IUNITR,  IS  USED. 
3  CONTINUE 

THE  LAST  READ  DID  NOT  ENCOUNTER  THE  END  CF  THE  FILE. 

READ  CIREA  0,1008  ) DUMMY 
1008  FORMAT  (A  1 C ) 


CHECK  TC  DETERMINE  IF  THE  PREVIOUS  READ  CF  UNIT,  IUNITR, 
ENCOUNTERED  THE  END  CF  THE  FILE.  SUBROUTINE  EOFILE  PERFORMS 
THIS  TASK, 
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SU3P.0UT  IFF.  RWOCRK  (IUNIT  E  *  N  PTSU  ♦  RH  0  ,  NPTSC  ,  RH  OC  ) 


THIS  $1)3  ROUT  I  HE  DETERMINES  IF  THE  RHO  ARRAY  IS  CORRECT. 

IF  IT  IS  NOT  THE  USc?  1$  SC  INFORMED  AND  THE  PROGRAM  IS  STOPPED. 


THIS  SUSP  CUTINc.  WAS  uAST  CHANGED  ON  *5  NOVEMBER  1981. 


DIMENS  IGF  RHO(NPTSU) ,RMOC (NPTSC) 


OATA  SMAlL/1. £-3/ 


1  CONTINUE 

FOR  THE  SAKE  OF  EFFICIENCY  SOME  LOCAL  V ARI A9LES  WILL  SE  SET  EQUAL 
TO  SOME  CF  THE  PARAMETERS. 

IUNIT  =  UNITE 
NPTSUL=NPTSU 
FiPTSCL=NF  TSC 

NPTSU  SHCULO  BL  POSITIVE  AND  EQUAL  to  NPTSC. 

IF(NFTSUl  ,LT.  0  .OR.  NPTCUL  . N E *  NPTSCL)  GO  TC  A 

2  CONTINUE 

DO  3  1  =  1, NPI SU L ,1 
RH  OI=RH  0  (I) 

RHOCIsRM)C(I) 

RCHECK= A  E  S (RHOI -RHOCI ) 

IF{ RCHEC K  .GT.  SMALL)  GO  TO  4 

3  CONTINUE 

ALL  THE  RHO  VALLES  AGREE  WITH  THE  CORRECT  VALUES, 

GO  TO  7 
*  CONTINUE 

WRITE  tlUMT,  10  0  4)  NPTSUL  ,  N  PTSC  L 
1034  FORMAT  (iH, 

11H  ,  2.6  HT  4  IS  IS  SUBROUTINE  RHOCHK.,/, 

21H  ,29HAF  ERROR  HAS  BEEN  DISCOVERED.,/, 

31H  , 57HTF  E  VARIABLES  WHICH  DEFINE  THE  RHC  ARRAY  ARE  NOT  CORRECT. 

4  ,2(1H  /), 

51H  ,  9H  FPTSU  =  ,  110,11  K,  OH  NFT3C  =  *I10,2tlH  ,/), 

61H  ,  3H  I  ,1Y,4>',  3HRHO,  5X,  11X,  AX,4HRM0C) 

5  CONTINUE 

N LIMITS  NPTS'J  L 
IF( NLIMIT  .LT,  NPTSCL)  NLIMIT= NPTSC L 
IF( NLIMIT  .lT.  1)  NLIMIT=1 

DO  6  N=  1 ,  NL  IMI T  ,  1 
WRITE!  IUN  17,10  0  0)  N,  R-M){  N  )  ,  RHOC  ( M) 

ICO  5  FORMAT (IF  , I  3, IX, L12.S, 1  1* ,El2. 5) 

6  CONTINUE 

WRITE ( IUN  IT,  1006) 

1005  FORMAT ( 1 F  ,/, 

liH  ,  43  HT  F  IS  SUBROUTINE  WILL  NEXT  STO3  THF  PROGRAM.) 

STOP 

7  CONTINUE 


ooo  OOO  OOOO  OOOO  OO  OOoOOOO  oooooo 


SUSROUTIf'E  FWC0N°  (IR  ■MT  APEN,CCMP,  NDIM,  MA<  PTE  »  MM3,  NP~U, 
1  WAVE.-  H) 


THIS  SUBROUTINE  R  ;  ADS  OR  WRITES  THE  AZIMUTHAL  COh  FF ICIENTS 
FOR  ALL  THE  P.HO  VALUES  AT  ONE  Z  LEVEL  FOR  ONE  FIELD  COMPONENT. 

IF  I RW  ECUALS  TWO  "HEN  A  WRITE  WILL  TAKE  PLACE,  CTHf  PH  IS"  A 
READ  IS  UOii.  EOF  THE  SAKE  OF  STYLE  AND  POSSIBLE  FUTURE  REVISION 
WHEN  A  READ  OPERATION  IS  DESIRED  CALL  THIS  SUBROUTINE  WITH  IRW 
EQUAL  TO  ONE. 


THIS  SUBROUTINE  WAS  LAST  CHANGED  ON  ?3  0CT03SR  19«1. 


COMPLc>  W AVEcH ( MDIM, MAXPTS ) 


1  CONTINUE 


FOR  THE  SAKE  OF  EFFIC IFNCY  SOME  LOCAL  VARIABLES  KILL  8E  SET  FOUAl 
TO  SOME  CF  THE  PARAMETERS. 

IRWL= IRK 
IT  AP£= IT  A  PEN 
MM3 L= MM 3 
NPT=NPTU 

NOW  A  TEST  WILL  BE  DONE  TO  DETERMINE  IF  A  READ  OF  A  WRITE 
HAS  BEEN  REQUESTED. 

IFfIRWL  .EG.  2)  GO  TO  12 

2  CONTINUE 

A  READ  IS  REQUESTED 

READ  (ITAFE, 1002  )  COM^ 

1002  FORMAT (29X, AID,/, 25X) 

3  CONTINUE 

READ  (ITAFE, 100  3  ) 

10C  3  FORMAT (76  X ) 

4  CONTINUE 

READ  (ITAFE,  100  4  )  (  <WA  VEE  H  ( I  ROW,  J  COL )  ,IROW  =  l  ,MN.3L,1>  »JCOL  =  i,NP"  ,  i) 
10Q4  FORMAT (5 (IX, E10. 4)) 

GO  TO  15 

12  CONTINUE 

A  WRITE  IS  REQUESTED 

WRITE(ITAPE,1012)  CC-HP 

1012  FORMAT ( 

11H  ,?8HTFE  NEAP  FIELD  COMPONENT  IS  ,Ai0,/, 

21H  ,  24HI (  FORMAT (2BX»A1C,/,25X)  ) 

13  CONTINUE 
WRIT£(ITAPE,1013) 

1013  FORMAT  ( 

11H  ,  48H  (  (WAVEEh  (IROW.JCOL)  ,IROW  =  l ,MM3, 1 )  ,JCOL=l  ,NPTU,1)  , 

2  27  H - IN  FORMAT  (5  (1X*E10.4)  )  ) 

14  CONTINUE 

WRITE (IT  A  PE, 10 14)  ( (WAVEE H ( IROW , JCOL ) , I ROW= \ , MM3L , 1) , JC  CL  - 1 , N°T , 1 ) 

1014  FORMAT (E ( 1H  ,E10.4) ) 
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SUBROUTINE  RWHEAC(IRN,ITAPEN,IUNITZ  . 

1  mA;;PTS,FF.EQMW,PCENT,  7INCH,  AOFI,  AOF*, 

2  MM3, NPT  U,RHO) 


THIS  SUBROUTINE  READS  OR  WRITES  THE  HEADER  TO  T  A  P  F ,  ITAPEN» 
IF  IRW  r*  DUALS  T  WG  TH_N  A  WRITE  WILL  7  A<  E  3LAC£,  OTHERWISE 
A  READ  IS  MADE.  FOR  THF  SAKE  OF  STYLE  ANC  FOSSIBLF  FUTURE 
REVISIONS,  WHEN  A  READ  OPERATION  IS  DESIRED  A  CALL  THIS 
SUBROUTINE  WITH  IRW  _OJAL  T0  ONE. 


THIS  SUBROUTINE  WAS  LAST  CHANGED  ON  ?7  DC  TO  3:  R  iBP-i. 


DIMENSION  RHOiNAXFTS) 


1  CONTINUE 

FOR  THE  SAKE  OF  EFFICIENCY  SOME  LOCAL  VARIABLES  WILL  BE  ;  ET  ECU. 
TO  SOME  CF  THE  PARAMETERS. 

IRWL= IRW 
ITAP£= IT  A  PEN 
IUNITsIUMTE 
MPTS=MA)F7S 
FMHZ=FSEQMH 
PC=  PCE  NT 
Z 1=  Z  IN  OH 
ai=aof: 

A  P=  A  0  F  P 
MM3  L=  MM3 
NPT=NPTU 

DETERMINE  IF  A  READ  OR  WRITE  OPERATION  IS  REQUESTED. 

IFCIRWL  .EQ.  2)  GC  TO  12 

2  CONTINUE 

A  READ  OPERATION  IS  REQUESTED. 

READ (IT AFC, 1002)  FMHZ,PC, ZI, AI, AP, MM3L , NPT 
100  2  FORMAT  (3eX,/,i5  (iX.Elt.4)  ,2  (IX.IIO)  ,/,i5X) 

A  TEST  WILL  3F.  PACE  TO  BE  SUPE  THAf  "'is  ARRAY,  RFG,  KILL  NCT 
BE  CLOBBERED  3Y  THE  NEXT  READ. 

I F ( M FT S  .LT.  NPT)  GO  TO  5 

3  CONTINUE 

READ  <ITAFt,iu33) 

1003  FORMAT ( 4 7 X ) 

-  CONTINUE 

READ (ITAFE- ' JG4: ’RHC  I),}  ,.NP7,i) 

1004  FORMAT (4  (1X,E10.4)> 

BEFORE  TFIS  SUBROUTINE  RETURNS  THE  f  IRA  METERS  FPEQMH,  FCENT, 
7INCH,  A  C  F I ,  AOFP ,  M‘!3,  AND  NPTU  KILL  BE  SET  EQUAL  TO  ?Ht 
VALUES  READ  FOR  FM-1Z.  PC,  ZI,  AI»  A  P,  **M3L»  AND  )  PT  CESPECTIVE  LY 


FR£0MH  =  F*'H7 


o  o  o 


PCENTr^C 
ZINCH=  Z I 
AOFI=A I 
ADFP=Ar 
MM3  =  Mt«3L 
NPTU=NPT 
GO  TO  1*5 
"  CONTINUE 

AN  ERROR  HAS  3EEN  DISCOVERED. 


WRITE(IUM7,1CC5>  IRWL,  F  MHZ ,  PC,  71 , A  I,  a  P  ,  MW3L  ,  NETS ,  NP~ 

1035  FORMAT  (1H, 2  (1H  ,/)» 

UH  »  26  l*TF  TS  IS  fUBRO  JT INE  RWHE4C,  »/ , 

Z1H  »2<JHAf  ERROR  HAS  BEEN  DISCOVERED.,/, 

31H  ,  55  NT  HE  NUMBER  CF  FHO  POINTS  IS  LARGER  THAN  THE  Ar'R  AY* 
4  2 (1 H  ,/), 

51H  ,9H  IRWL  =  ,115  ,/. 

61H  ,  9HFW  E  CNH  =  ,E  1  5.  3 , 5.V ,  9H  PCE NT  =  ,E15.9,/, 

71H  ,9H  ZINCH  =  , El 5.3:,/, 

S1H  ,  9H  A  CFI  =  ,£ 15.9 , 5X  *  9H  AOFP  =  ,El5.9,/, 

91H  ,9H  MM3  =  ,115  ,5X»9HMAXFTS  =  ,115  ,/, 

01H  »9H  MPT  =  ,115  , 2 ( 1 H  ,/) , 

11H  ,  42  H”  (-  IS  SUBROUTINE  WILL  NOW  STOP  THE  ’P.OGRaH.) 

STOP 

12  CONTINUE 

A  WRITE  OPERATION  IS  REQUESTED. 

WRITE  ( IT  A  Pt,  10  12)  FMHZ,PC.ZI,AI,AP,MM3L.N*»T 

1012  FORMAT  ( 

11H  ,  37H.FR  EQMH,  PCFNT ,  ZINCH ,  AOF I,  AOFP  ,  MM3  ,NP7U  , /, 

2  5(1H  ,E10.4),2(H  ,110),/, 

31H  , 44HI f  FORMAT(3«X,/,5 (iy,E10.4) ,2(1Y,I10) ,/,45X)  ) 

13  CONTINUE 
WRITE(ITA°E  ,1013) 

1013  FORMAT  ( 

UH  ,  46H  (RHO(I)  •  1=1,  N97U  ,  1  >  —  IN  FORMAT  <  4  (IX  *E10  ,4  i  )  ) 

14  CONTINUE 

WRITE  { IT  A  PE,  10 14)  (P'-iOd  )  ,1  =  1,  N°T  ,1 ) 

1014  FORMAT  (4  (li  ,ElC.  4)  ) 

15  CONTINU.-. 

RETURN 
END 


O  O  o  o  o  o  o  O  oo  no  OO  OOOOO  OO  OOOOOoOOOOOOOO 


SUSRO'JT  IN  E  RW1W1Z  CRW, ITAPEN,  IUNITE* 

1  MAXPTS ,FREQMH,PCEN",  Z INCH, A  OF  I, AOF 3, 

2  MM  T  »  NPTU »  P.HO  »  MDIM  ,  'IAV  EE-U 


THIS  SO’FOUTINF  ALONG  WITH  ITS  TREE  ?EADS  °p-  WMTES  THE 
INFORMATION  ON  TA  PE,  IT  A  PEN*  CONCERNING  THE  INCIDENT  “LANE 
WAVE  AND  THE  AZIMUTHAL  COEFFICIENTS  FOR  ALL  ~HL  RHO  VALUES 
AT  THE  Z  ZLEVFL,  ZINCH.  IF  IRW  EQUALS  TWO  THEN  A  WFITE  WILL 
TA<E  PLACE,  OTHERWISE  A  READ  IS  MADE.  FOR  THE  SAKE  OF  ST  >  LE 
AND  POSSIBLE  FUTURE  REVISION  WHEN  A  READ  OPERATION  IS  DESIRED 
CALL  THIS  SUBROUTINE  WITH  IRW  EQUAL  TO  ONE. 


THIS  SUBROUTINE  WAS  LAST  CHAFED  ON  25  OCTOBER  1981. 


COMPLEX  WAVEEH{MDIM, MAXPTS  ,6) 

DIMENSION  RHO(  MAXPTS  ♦  ,C0MPPN<6)  ,C0MPWN<6> 


THE  ELEMENTS  OF  THE  ARRAY  COHPWN  STORE  THE  NAME  CF  EACH  CIELD 
COMPONENT . 

DATA  C  CM  F WN/lOhFPHI  , 10HERHO  ,10HEZEE  , 

1  10HHFHI  , 1CHHRHO  ,13HHZEE  / 


i  CONTINUE 

FOR  THE  SAKE  OF  EFFICIENCY  SOME  LOCAL  VARIABLES  W IL  L  BE  SET  EQUAL 
TO  SOME  CF  THE  PARAMETERS. 


IRW L= IRW 
IT  APE= IT  A  PEN 
XUNITs  2UMTE 
MPTS=  MA>  °TS 
FMHaFREOKH 
PC=  PCS  NT 
ZINL=  ZINCh 
AI= AOF  l 
AP=AOFP 
MM3 L= MM  3 
NPTU  L=  NPTU 
MD IML=  MO I M 

2  CONTINUE 

THE  HEADING  WILL  NCW  Br  READ  OR  WRITT  EN .  3U3P01TINE  RWFE  AO 
WILL  PERFORM  THIS  TASK. 

CALL  RWHEAO(IRWL,ITA’F.,:UNIT, 

1  MPTS, FMH.PC, Z1NL, AI.AP, 

2  HM3L,  N^T  jl,  RHO) 

A  TEST  WILL  3E  CuNE  TO  BE  SUCE  THAT  THE  FIRST  DIMENSION  CF 
THE  ARRAN,  WAV  EE F,  IS  LARGE  ENOUGH. 

IF(MM3L  .GT.  3  .ANT.  HMSL  .LE.  MD1MU  GO  TO  i 

3  CONTINUE 


AN  ERROR  HAS  3EEM  DISCOVERED,  SC  THE  USER  WILL  BE  INFORMED 
AND  THE  F  ROGRA.M  STOPPED* 


WRITE  ( IUMT*iO  C  3)  IPNL. ITA  DE , MPTS,  NPTU  L ,  F 1H ,  PC*  A I ,  A  P, Z INL  ♦ 

1  MM3L.  MO  I  ML 

1903  FORMAT  (1H, 2(1H  ,/), 

11H  ,26hTHS  IS  SUBROUTINE  PW1W17.,/. 

21H  , 52  H  T  r  E  NUMBER  OF  AZIMUTHAL  COEFFICIENTS  FOR  ~HE  CURRFNT  , 
31H  »  52 HZ  LEVEL  SHOULD  BE  POSITIVE  BUT  LFSS  THAN  CR  EOUAL  TO,/. 
41H  ,  53  HTF  £  FIRST  DIMENSION  CF  THE  A  R°A  Y ,  WAVEFF.  THIS  IS  N  CT ,  / 
51H  ,  3 HI K:;  CASE.,2(1H  ,/), 

61H  ,  9H  IRW  =  ,110  9HITAPEN  =  ,110  ,/, 

71H  ,  9HMAXPTS  =  ,110  ,11X,  9H  NPTU  =  ,110  ,/v 

91H  ,  9HFFEQMH  =  ,  E 10 . 4, 1  lX  ,  9H  PC ENT  =  ,El0.4,/, 

91H  ,  9H  AOFI  =  «  E10.  A,  1 IX,  9H  AOFP  =  ,E10.W, 

OH  ,  3 H  ZINCH  s  .L10.4,  /, 

llH  ♦  9H  MM3  =  ,110  , 1 IX ,  9H  MOIM  =  ,110  ,2(1H  ,/), 

21H  ,  42  RTF  IS  SUBROUTINE  WILL  NON  STOP  THE  PROGRAM,) 

STOP 

4  CONTINUE 

A  TEST  WILL  BE  DONE  TO  DETERMINE  IF  a  READ  0*  A  fcRITF 
HAS  BEEN  REQUESTED. 


IF  ( I FWL  .EO,  2)  GO  TO  15 
5  CONTINUE 

A  READ  OPERATION  HAS  3EEN  REQUESTED. 

NOW  THE  AZIMUTHAL  COEFFICIENTS  F CR  EACH  COMPONENT  OP  THE 
NEAR  FIELD  WILL  3E  READ.  THE  REPEATED  CALLS  TO  SI  BP OUT  INE 
RWCOMP  PERFORMS  THIS  TASK. 

DO  6  I  CO t °=1, 6,1 

CALL  RWCC  MP  (IRWL,  ITAOE,  COMPRNt  ICOMP)  ,Mf)  XML  ,  MPTS  ,  MM3L ,  NFTUL , 

1  WAVEEHtl.l, ICOMP)) 

o  CONTINUE 

THE  COMF1 NENT S  SHOULD  HAVE  BEEN  READ  IN  THE  ORCFR  THAT  THEY 
ARE  LISTED  IN  The  ARRAY ,  COMPWN. 

DO  7  ICON P=i, 6,1 

IF  (CCMP'-N  (ICOMP)  ,NE .  C  OMPWN  (ICCMP)  )  GO  TO  3 
7  CONTINUE 

THE  PARAMETERS  FREOM«tPCENT,ZIUCH,ACPI,AOFP,MMS,AND  NPTU  WILL 
SET  EQUAL  TO  THE  VALUES  READ  FOR  M  PTS,  p  MH ,  PC .  Z  J.NL  ,  A I ,  A  P ,  mM  3L , 
NPTU  L  RESPECTIVELY. 

FREQMH=FMH 
POE  NT  =  PC 
ZINCH=ZINL 
A0FI= A  I 
AOFp=AF 
MM3  =  MM  3L 
NPTJ=NF  TLL 
GO  TO  l7 
)  CONTINUE 

AN  ERROR  PAS  3FFN  L ' ..COVERED.  THE  CCMPCNFNTS  WERF  NOT  RF AC 
IN  THE  CRGfcF.  EXPECTS).  THE  USER  WIlL  9E  SO  IMFCRMEC  AND  ThEN 


o  o  o 


C  THE  PROGRAM  STOPPED. 

C 

WRITE  CUMT.iOOc)  :R*L.  ITAPE,MCTS,  FFTUl,FMH,  PC,  A  I , A  P,  1  INI » 
1  MM3 L*  MO  I  ML 

1003  FORMAT  (1H,2(1H  ,  /)  , 

HH  »  26 FTP  IS  IS  SUBROUTINE  RW1W1Z. ,/, 

21H  ,  51  HT  P  C  COMPONENTS  WERE  NOT  READ  IN  THE  EXPECTED  CRCE«. 
3  2(iH  ♦  /), 

41H  ,  9H  IRW  =  ,  I1C  ,11X,  9HITAPEN  =  ,110  ,/ , 

51 H  ,  9HMAXPTS  =  .110  ,11X,  9H  NPTU  =  ,110 

61H  ,  9HFFEOMH  =  ,Ei3.4,ilX,  9H  PCENT  =  ,E13.A,/, 

71H  ,  3H  AOFI  =  ,E1G.4,11X,  9H  AOFP  =  ,E10.4,/, 

91H  ,  9P  7  INCH  =  ,E.10.4,  /, 

01H  ,  9H  MM3  =  ,I1C  ,  1 IX ♦  9H  MDIM  =  ,110  ,2<1H  ,/), 

11H  ,  3H  I  « IX, 14H  EXPECTED  ORDER , IX, 1CHORDER  READ) 

9  CONTINUE 

00  10  1=1, 6,1 

WRITE  ( IUMT « 13  DQ)  I,COMPWN(I)  ,CUMPRN  (I) 

1009  FORMAT  (1 P  , IX , 1 1, IX , IX, 2X , A  10 * 2X, IX , A 1 0 ) 

13  CONTINUE 

WRITE(  IUMT ,1310) 
iOlQ  FORMAT  (2  <1H  ,/)  , 

11H  , 42  FTP  IS  SUBROUTINE  WILL  NOW  STOP  THE  =»ROGPAN.) 

STOP 

15  CONTINUE 

A  WRITE  OPERATION  WAS  REQUESTED. 

DO  16  ICC MP=1, 6 ,1 

CALL  RWCC  f'P  ( IRWL,  IT  A  3E»  C  OMPWN  (  ICOMP )  ,M3  IFL,MPTS,  1M31  ,  NP’UL 
1  WAV EE P (1,1, ICOMP)  ) 

16  CO  NT  IN  IE 

17  CONTINUE 
RETURN 
EN3 


OOO  OOO  OOOOOOOCTOOOOOOOOO 


SUBROUTINE  TCHECK  (IL'NITE  ,FMHZW.FMHZG,  pCFNTW.pCENTG) 


THIS  SUBROUTINE  COMPARES  FMHZW  TO  FNHZG,  AND  FCENTW  TO  FCENTG « 
IF  BOTH  COMPARISONS  -AiRE  CLOSE  THEN  THIS  SUBROUTINE  PETUr.NS  TO 
THE  CALLING  SUBPROGRAM.  OTHERWISE  IT  SO  INFORMS  THE  USER  AND 
THEN  STOPS  THE  PROGRAM, 


THIS  SUBROUTINE  WAS  LAST  CHANGED  ON  5  NOVEMBER  1981, 


DATA  SMALL/1, E-3/ 


1  CONTINUE 

FOR  THE  SAKE  OF  E  FFICI£MCY  SOME  LOCAL  VARIABLES  WILL  PE  SET  EQUAL 
TO  SOME  CF  THE  PARAMETERS. 

IUNIT=IUMTE 
FW=FMFZW 
FGsFMFZG 
PC W=  PCE  NTW 
PCG=  FCENTG 

2  CONTINUE 

NOW  THE  COMPARISIONS  WILL  BE  MADE. 

A3SFD  =  AES(FW-FG) 

IF  C  ABSFC  *GT.  SMALL)  GO  TO  5 

3  CONTINUE 
ABSPCD=AES(PCW-PCG) 

IF( A  PS  PC  C  .GT,  SMALL)  GO  TO  5 

4  CONTINUE 
GO  TO  6 

5  CONTINUE 

WRITE (IUMTf  10  0  5)  FW,  PCW,  FG,PCG 
1005  FORMAT  (1H,/, 

11H  « 26HTFTS  IS  SUBROUTINE  TCHECK.  ,/, 

21H  ,29HAN  ERROR  HAS  BEEN  DISCOVERED.,/, 

31H  ,4/HTHE  FREQUENCY  OF  THE  INCIDENT  PLANE  WAVE  CP.  THE,/, 

41H  , 49  HpE  RCENT  AGE  OF  WATER  CONTENT  IN  THE  SOIL  READ  FROM,/, 

51H  , 51 HTFE  FILE  WITH  THE  AZIMUTHAL  COEFFICIENTS  IS  NOT  THE,/, 

61H  » 45  HS  A  ME  AS  WHAT  THE  USER  SPECIFIED  IN  THF  INPUT, ,2(1H  ,/), 

71H  ,17X,2X,2<, 12H  FREQUENCY  IN, ?X,?X ,lX, 13HPSRCL NTAGE  OF,/, 

81H  »17X»2X»3X* i OH  MEGA  HERTZ, 3X, 2K , IX, 1 3HWA TER  CONTENT,/, 

91H  , 17HS  FECIFIEC  BY  USER ,  2X  , IX ,  E14.  8,  IX  ,  2X , IX  ,E  14  .8  , /, 

OlH  ,17H  IN  TH>  FIL.  ,2X,1X,E14,  3,lX,2<tlX,El4.8,Z<  IH  ♦/), 

11H  ,J.3hTHS  SUBROUTINE  WILL  NEXT  STOP  THE  PROGRAM.) 

STOP 

5  CONTINUE 

THE  FRLOLcNCY  AND  TH-_  PERCENTAGE  ARE  AS  DESIRED. 

RETURN 

END 
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OOOO  OO  OO  OOOOOOOOOOOOOOOOOOOOO 


SU3R0UTIHE  WAZICt  IL'NTTO,Hr  IM.RHO, IFLAG. MUSED, 
i  M.AGFI.AOFP, ZINCH,  C0MPN1.C0MPN2, AZTC) 


THIS  SUBS  OUT  INF  WRIT  .S  TO  UNIT,  IUNITO,  FOR  ALL  Thf  FHO  VALUES  • 
ONE  CF  THE  USER  SPECIFIED  SETS  OF  AZIMUTHAL  COE FF ICTENTS .  IF 
IFLAG  EQUALS  ONE  HE  COEFFICIENTS  ARE  OF  ONE  01-  THE  COMPONENTS 
OF  THE  SCATTERED  ELECTRIC  FIELD.  IF  IFLAG  OOi'S  NCT  t  DUAL  ONE  THE  H 
THE  COEFFICIENTS  ARE  OF  THE  PRODUCT  OF  ETA  0  AND  ONE  CF  THE 
COMPONENTS  OF  THc  SCATTERED  MAGNETIC  FIELD. 

THE  CALLING  SUBPROGRAM  IS  RESPONSIBLE  FOR  A  CONSISTENT  SETTING  CF 
THE  FAR  A  f  ETERS  IFLAG.  COMPNi,  AND  COMPN2. 

NOTE  NO  ERROR  CHECKING  IS  PERFCRMFO. 

THE  OUTPIT  TO  UNIT.  IUNITO.  WRITTEN  BY  THIS  SUBROUTINE 
IS  SUITAELE  TO  EE  PRINTED  ON  A  LINE  PRINTER. 


THIS  SUBROUTINE  WAS  LAST  CHANGED  CN  14  JANUARY  1982. 


COMPLEX  AZICIMOIM.  16) 


DIME  NS  10 1'  RH0(16>.  U  MUSED) 


1  CONTINUE 

FOR  THE  JAKE  OF  EFFICIENCY  SOMz  LOCAL  VARIABLES  WILL  BE  SET  EQUAL 
TO  SOME  CF  THE  PARAMETERS. 

I'JNIT  =  IINIT0 
IFLAGL=IFLAG 
MTOT  AL=  Ml SEO 
CNAME1  =  CCMFN1 
CNA  ME2  =  CC  NPN2 
C 

IF ( IFL AGL  .Nil.  i)  GO  TO  3 
C 

2  CONTINUE 
C 

KRITEdUMT,  1002)  CNA  HE  1 

1002  FORMAT  (iH./.lH  .57X.A5) 

C 

C-0  TO  A 
C 

3  CONTINUE 
C 

WRITE(  IUMT.1J03)  CNA  ME  1 

1003  FORMAT  (l  H./.1H  ,54X,AlO> 

C 

*  CONTINUE 
C 

WRITE  ( lUMT.iOOD  AOFI,  AO FP,  ZINCH 

1Q3  4  FORMAT  < 

11H  »  30X. 5 1H7HE  ANGLE  OF  INCIDENCE  OF  THE  PLANE  WAVc  IN  DEGREES, 

2  -.H  IS  ,F4. 1 ,  /, 

31H  »  30  X  ♦  S  IhT  H  E  POLARI7ATION  AN  GLr  OF  THE  PLANE  WAVE  IN  DEGREES  » 


4  4H  IS  , F4.  i  «  /, 

51H  , 21X  »  £ 1HTHE  SCATTERED  FIELD  IS  EVALUATED  ABOVE  THE  GPCLND  A, 

6  23H  DISTANCE  IN  INCHES  OF  ,F3*1*2(1H  ,/)) 

C 

IF( IFLASL  .ME.  1 )  GO  TO  6 
C 

5  CONTINUE 
C 

WRITE  (IUN  IT,  10  0  5)  CN4ME  2  »C NA  M£? ,  CNA  MT'2 , CNA  MS? 

IOC 5  FORMAT ( 

11H  »  33  FTH  E  AZIMUTHAL  COEFFICIENTS  OF  TH  E  ,  A  6 , 1 2r  CCMPONEM  OF, 

2  43 H  THE  SCATTERED  FIELD  ARE  LISTED  UNDER  THEIR, 

3  23H  CORRESPONDING  VALUE  OF,/, 

41H  , 46HRH  0  IN  UNITS  OF  VOLTS  PER  METER.  TO  OBTAIN  THE, A6, 

5  3SHCC  MPONENT  OF  THE  TCTAl  FIELD  ADD  THE.A6,  12HCGMPCNENT  CF, 

6  13H  THE  INCIDENT,/, 

71H  , 1 3HF I ELO  AND  THE , A5, 35H00Mo0NENT  OF  THE  FItLt  SCATTERED  aY, 

8  50 H  THE  AIR  GROUND  INTERFACE  TO  THE  SUM  OF  THE  SFRIES, 

9  11 H  CEFINED  BY,/. 

OlH  ,  56H7HESE  AZIMUTHAL  COEFFICIENTS.  THE  MAGNITUDE  WILL  TRFN  BE, 

1  57H  F  OR  MALI  ZED  FOR  AN  INCIDENT  PLANE  WAVE  HAVING  AN  ELECTRIC- 

2  6H  FIELD,/, 

31H  ,22H0F  ONE  VOLT  PER  METER.) 

C 

GO  TO  7 
C 

6  CONTINUE 
C 

WRITEdUMT.lflCe)  CNAME2  ,  CNA ME2  ,  CNA  FE2 ,  CNAME2 
10 3  o  FORMAT ( 

11H  ,  33  HTF  h  AZIMUTHAL  COEFFICIENTS  OF  THE  °  FODUC  T  OF  ETAo  (IN, 

2  14H  CHMS)  AND  THE,  A  6 , 3  6HC0M  PONENT  CF  THE  SCATTERED  FIELD  ARE, 

3  7 H  LISTED,/, 

4lH  ,  56 HU  f  DER  THEIR  CORRESPONDING  VALUE  OF  RHO  IN  UMTS  Cc  VOLTS, 

5  25 H  FFR  METER,  TO  OBTAIN  THE, A6 , 22HCOMPONENT  OF  THE  TOTAL, 

6  10H  FIELD  ADD,/, 

71H  ,  3HTFE,A6, 3  9H  COMMON  ENT  OF  THE  INCIDENT  FIELD  AND  TFE.A6. 

6  53  HCC  MPONENT  OF  THE  FIELD  SCATTERED  BY  THE  AIR  GROUND, 

9  13H  INTERFACE  TC,/, 

OlH  , 51 HTHE  RATIO  OF  THE  SUM  OF  THE  SERIES  DEFINED  BY  THESE, 

1  57H  AZIMUTHAL  COEFFICIENTS  AND  E^AO  (IN  CHHS) ,  THE  MAGNITUDE, 

2  10  H  WILL  THEN,/, 

3lH  *  50  HBE  NORMALIZED  FOR  AN  INCIDENT  PLANE  WAVE  FA  V  l'NG  AN, 

4  38 H  ELECTRIC  FIELD  OF  ONE  VOLT  PER  METER. ) 

C 

7  CONTINUE 
C 

DO  10  1  =  1,  u,l 
II-(4* (1-1) ) fl 
II»laII*l 
IIP2=IIf2 
IIP3=I If 3 

WRITE  ( IUMT,10  0  7)  RHO  (II)  ,  RHO  ( I  IP  1)  ,RHC  (IIP2)  ,RHO  ( 1 1 P 3  ) 

1007  FORMAT  (2  (1H  ,/, ), 

11H  »  7X » 9HAZIML‘TH A\_,  *» (2X  , 6FRH0  =  ,F4.1,12H  CENTIMETERS),/, 

21H  ,  7X , 9  H  ORDr  P  , 4 (2X , 3X , 4MREA L , 3X , IX, 10H IMA  GINA R v  ,lX)) 

C 

■)  CONTINUE 
C 

DO  9  IM=  1  ,NTDT  A  L,  1 

WRITE(  II  MT,  10  CC  )  M(  IM)  ,  AZIC  ( IM  ,11  )  ,  A  Z IC  (  IN ,  I  IP  1 )  ,A7IC  ( IM  ,  I  IP?)  , 

1  AZIC (IM,II?3) 
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Appendix  III 


Table  of  Symbols 

v-, 


O 


0 


tl; 


C 


I1  Electric  field  vector  of  the  incident  plane  wave 

El  j  component  of  the  electric  field  of  the  incident 

J  plane  wave,  where  j  =  x,  y  or  z 

H1  Magnetic  field  vector  of  the  incident  plane  wave 

Hi  j  component  of  the  magnetic  field  of  the  incident 

J  plane  wave,  where  j  =  x,  y  or  z 

"Ep  Electric  field  vector  of  the  primary  field 

Ep  e;  component  of  the  primary  electric  field,  where 
Q  t,  =  x,  y,  z,  p  or  <j> 

Hp  Magnetic  field  vector  of  the  primary  field 

Hp  c  component  of  the  primary  magnetic  field,  where 
;  C  =  x,  y,  z,  p  or  <j) 

Es  Scattered  electric  field  vector 

E^  c  component  of  the  scattered  electric  field,  where 
Q  K  =  x,  y,  z,  p  or  <p 

H5  Scattered  magnetic  field  vector 

H*  c  component  of  the  scattered  magnetic  field,  where 
^  C  =  x,  y,  z,  p  or  <j> 

"E  total  electric  field  vector 

E  c  component  of  the  total  electric  field 

Q  C  -  x,  y,  z,  p  or  <j> 

H  total  magnetic  field  vector 

H  x,  component  of  the  total  magnetic  field 

5  C  =  x.  y,  z,  p  or  $ 

<J>1  First  Azimuthal  potential 

i^2  Second  Azimuthal  potential 

e  base  of  natural  logarithm 


1 
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£q  permittivity  of  free  space 

eg  Relative  dielectric  constant  of  earth 

em  Relative  dielectric  constant  of  land  mine 
e.j  imaginary  part  of  a  relative  dielectric  constant 
real  part  of  a  relative  dielectric  constant 
tan  6  Loss  tangent  of  a  dielectric  material 
a  attenuation  constant  of  a  dielectric  medium 

kQ  wave  number  of  free  space 

£  propagation  vector 

n0,ETAO  intrinsic  impedance  of  free  space 
\Q  free  space  wave  length 

f  frequency 

c  speed  of  light 

p  reflection  coefficient,  parameter  of  cylinderical  coordinate 
system 

analogue  for  earth  of  transmission  line  load  impedance 
ZQ  analogue  for  earth  of  intrinsic  transmission  line  impedance 
Re{*}  real  part  of  the  bracketed  expression 
Im( * }  imaginary  part  of  the  bracketed  expression 


